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FROM GAP PROBABILITIES IN RANDOM MATRIX THEORY TO EIGENVALUE 

EXPANSIONS 

THOMAS BOTHNER 


Dedicated to Percy Deift and Craig Tracy on the occasion of their 70th birthdays 


Abstract. We present a method to derive asymptotics of eigenvalues for trace-class integral operators 
K : L^(J; dA) 0, acting on a single interval J C M, which belong to the ring of integrable operators [28]. 
Our emphasis lies on the behavior of the spectrum {Ai(J)}^Q of X as | J| ^ oo and i is fixed. We show 

that this behavior is intimately linked to the analysis of the Fredholm determinant det(7 — as 

|J| ^ oo and 7 t 1 a Stokes type scaling regime. Concrete asymptotic formulae are obtained for the 

eigenvalues of Airy and Bessel kernels in random matrix theory. 


1. Introduction and statement of results 

A well-known result in orthogonal polynomial random matrix ensembles states that the spacing distribu¬ 
tions of eigenvalues oi N x N Hermitian matrices are encoded in Fredholm determinants [31, 23]. In more 
detail, for a given positive weight function w{t)^ the probability E]s[{n; J) that a matrix from the ensemble 
associated with w{t) has n C Z>o eigenvalues in J C M equals 

EAr(n; J) = det(/-7AAr)|^^^. (1.1) 

Here, is the finite rank integral operator on L^(J;dA) with kernel 


iV-l 

^ Pj e C[t] : 

j=o 

The Fredholm representation of eigenvalue spacing statistics is also valid in the limit A —>■ oo and the 
resulting trace-class integral operator K : L^(J;dA) O depends crucially on the Hermitian model we start 
with and at which local point in the spectrum the scaling is considered. 

For instance, two of the most commonly encountered kernels in random matrix theory arise in the Gaussian 
Unitary Ensemble (GUE) by scaling in the “bulk” of the spectrum [23, 34, 39], resp. at the “soft edge” 
[10, 26, 32]: 


Pj{t)pk(t)wit)At = Sjk- 


Asin(A, p) 


KAi{X,p) 


sin(A — p) 
tt{X-p) ’ 
Ai(A)Ai'(/r) 




(-s,s). 


Ai'(A)Ai(/i) 


X — p 


s > 0; resp. 


Jai = (s, oo), s G K; 


( 1 . 2 ) 

(1.3) 
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i.e. the sine kernel and the Airy kernel, where Ai(z) is the Airy function. A third well-known kernel appears 
when one scales the Laguerre or Jacobi Unitary Ensembles (LUE or JUE) at the “hard edge” [26, 39]: 






Asess (b; S ^ 0, CL ^ 1, 


(1.4) 


2(A-m) 

i.e. the Bessel kernel, where Ja{z) is the Bessel function of order a. Erom the random matrix point of view 
one is interested in the analogue of ( 1 . 1 ), 


C_l\n jn 


7=1 > 


L>(J;7) = det(/-7A:)|^2(j) 


which equals the probability that there are n G Z>o bulk scaled {K = K^in), resp. soft-edge scaled {K = K^i), 
resp. hard-edge scaled {K = ATg^gg) eigenvalues in the interval Jsm, resp. Jai, resp. Jeess-^ 


We shall denote the eigenvalues of the integral operator K by {Ai(J)}“Q. In case of the three kernels 
(1.2), (1.3) and (1.4) it has been proven [35, 27, 36, 37] that the spectrum is simple and we can order 

1 > Ao(J) > Ai(J) > ... > 0. (1.5) 

The lower bound follows from positivity and the upper one from the fact that K : L^(IR;dA) O with 
K = Arsin,ArAi and K : L^((0,oo);dA) O with K = are projection operators (cf. [35, 36, 37]). In 

particular, the last statement implies by minimax characterizations that for each fixed i G Z>o we have 

A,(J)^1, as [Jl^-boo Ai(J)^l, as K = 

The eigenvalue behavior (1.5) provides us with a quick qualitative explanation of the phase transition near 
7 = 1 in the asymptotic behavior of D[J-,"f) as [J] —>■ oo, cf. [40]. Indeed, use Lidskii’s Theorem and write 

OO 

A>(J;7) = det(/-7A:)|^2^jj = (l - 7Ai(J)), (1.7) 

i=0 


SO that (1) for 7 < 1 all factors are bounded away from zero. Hence D{J',j) approaches zero exponentially 
fast. On the other hand (2) for 7=1 the i-th factor will approach zero for each i and D{J-,^) tends to 
zero faster than for 7 < 1. In the remaining case (3) for 7 > 1 the determinant will vanish for a discrete set 
{dijigz, the solutions of Xi{J) = 7 “^. 

As another application, the eigenvalue asymptotics (1.6) can be used to derive quantitative asymptotic 
information for the normalized eigenvalue probabilities 


r(u; J) 


A(u; J) 


E 


(1_A,J.....(1-A,J’ 


n G Z>i, 


i.e. (1.6) serves in the exact evaluation of “large gap probabilities”, see [36, 37] for concrete formul®. For 
this to work one requires subleading terms in ( 1 . 6 ) and a standard way to obtain these is the method of 
commuting differential operators: by a remarkable coincidence the three integral operators on A^( J; dA) with 
kernels (1.2) and (1.3), (1.4) commute with the operators given by 


Csinif] = 



f, 


and 


/’(a) 

'^Bess 


[/] 





( 1 . 8 ) 


defined on appropriate function spaces. Since (i^, C) share the same eigenfunctions the desired eigenvalue 
expansion is then obtained from WKB arguments applied to the differential equation. This is exactly the 


^Equivalently, £'(0; Jai) a-nd £(0; Jbg 
from the GUE and LUE. 


5 ) are the distribution functions of the (rescaled) largest and smallest eigenvalue drawn 
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approach employed by Poliak, Slepian and Fuchs [35, 27] for i^sin: the eigenfunctions of £sin are prolate 
spheroidal wave functions, cf. [33], and the following asymptotics was derived, 

l-A,(Jsin) = ^23*+"s*+Je-2«(l + 0(s-i)), s^+oo, (1.9) 

valid for fixed i € Z>o. 

In case of (1.8) no explicit formulm for eigenfunctions are known, still Tracy and Widom [36, 37] used 
the operators (1.8) in combination with a trick and derived analogues of (1.9), see Corollaries 1.7 and 1.10 
below. This derivation was somewhat heuristic and it seems to work well only for Sturm-Liouville operators 
£ with rather simple potentials. For more general kernels in random matrix theory, compare Section 1.3 for 
a short discussion, the potentials would be transcendental functions of Painleve type and the method fails. 

For this reason we choose a different approach to the asymptotics of Ai(J) as j Jj —>■ oo which will only 
rely on the integrable structure [28] of the kernels, i.e. them being of the form 

A-(A.rt = (Lio) 

A — /i 

The idea is to obtain “large gap asymptotics” for D{J;j) as j Jj —)■ oo and simultaneously 7 f 1 in a specific 
scaling regime. This regime has to be chosen in such a way that we can observe individual factors of the 
infinite product (1.7), i.e. preferably we would like to produce an asymptotic expansion of D(J;^) which 
involves a finite product with distinct factors. 


Remark 1.1. It is not clear at all how to read ojf individual factors from the expansions of D{J;"f) as 
jjj —>■ 00 and 'y < 1 is kept fixed: For instance, for the sine-kernel determinant [2, 12, 30, 24, 18], as 
s —)■ + 00 , 


£’(Jsin;7<l) = exp 


Fl( Jsini 1 ) — exp 


2v 

- s 

TT 


(4,)^G2(1 + ^)G^(1-£)(1 + 0(s-i)), u = -ln(l- 7 ); 


s ^Co(l + O(s ^)), 


2 ?! 

Co = exp 


-In2 + 3C'(-1) 


( 1 . 11 ) 


involving Barnes G-function G(-) and the Riemann zeta-function ((■)■ 


In order to motivate the usefulness of the phase transition point 7=1 we recall the following result from 
[ 6 ]. This result was obtained with the help of the known behavior (1.9), i.e. the line of reasoning is reversed 
in loc. cit. - nevertheless it will provide us with valuable insight for the general case ( 1 . 10 ) where little 
information on Xi{J) is available. The result we are referring to is the Theorem below. ^ 


Theorem 1.2 ([ 6 ], Theorem 1.12). Given x S K /et p = p(x) G ^>1 such that p = 1 for x < 5 
X +5 <P<X+i for X ^ 5 - There exist constants sq = so(x) > 0 ond vg = uo(x) > 0 such that 


■D(Jsin;7) = exp 


.-con(i 

2=0 


^—v '^^('^sin) 

1 ('^sin) 


(l + o(, 


-min{p-x-5.1} 


( 1 . 12 ) 


uniformly for s > so,v =— j) > vq and v > 2s — x^t^s. The universal constant cq appeared in (1.11). 


Note that from (1.9), we have explicit information on the factors, 

1 + e-” , = 1 + A2-3*-2s-*-5e2^-” (1 + 0 (s-^)) , 0 < * < p - 1 

l-Ai(Jsi„) y/TT 

and these contribute in general to the leading behavior of D{JsinTX) 9'S J+sinj —>■ 00 . In fact by construction 
of p G Z>i we have 0<p — x~ 5 <l and thus for w = 2s — x In s, 

^ = p-l: Ap_i(J,i„))-' ~ = o(I); 

z = 0,... ,p - 2 : e-"(l - A,(J,i„))-' ~ + o(I) 


^We have changed the definition oi v = —^ ln(l — 7 ) 
kernels considered in this paper. Also, [ 6 ] uses Agin (A,/i) 


in [ 6 ] to n = — ln(l — 7 ). This allows us to have a uniform v for all 
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SO that the asymptotics of D{Jsin',j) changes each time we cross one of the infinitely many Stokes curves 

u = 2s-xlns, x=P-^, peZ>i- 

In short, the asymptotic Stokes phenomenon of the Fredholm determinant encodes the information of the 
spectrum we are after. The natural idea is now to reverse the procedure which lead to (1.12), i.e. 

(A) Identify the asymptotic Stokes region for a given determinant and derive transition asymp¬ 

totics of type (1.12), i.e. derive an expansion of the form 

p-i 

+ |J|^+oo, 7tl, 

i=0 

with distinct factors fi{J) and where p S Z>i is tailored to the double scaling regime. 

(B) Extract from fi{J) the behavior of Ai(J) as |J| oo and i G Z>o is fixed. 

Remark 1.3. The existence of a phase transition for D{J\x) near 7 = 1 is a priori ensured for many, if 
not all, kernels in orthogonal polynomial random matrix theory. The operators are positive definite and since 
D{J; 1) has a probabilistic interpretation we have automatically 

1 > Ao(J) > Ai(J) > ... > 0. 

In this paper we carry out both steps for and discuss applicability of the results to a more 

general Painleve type kernel to be discussed below. 


1.1. Results for the Airy and Bessel kernel. The following two Theorems are the major results of the 
manuscript. 


Theorem 1.4. Given x G K determine p = p(x) G ^>0 such that p = 0 for x < “5 o.n-d X +5 <P<X+| 
for X > There exist positive to = to(x) ^ = '*^o(x) such that 


D{JAi]l) = exp 


p-i 


S 8 


® % n (^+ 




\f7^ 


(l + o(i 


-min{p-x-i 


■5>)) (1.13) 


uniformly for t = (—s) 2 > to,v = — ln(l — 7 ) > uq and v > %'J‘^t — xln^- Here 


rjo = exp 


-In2 + C'(-1) 


involving the Riemann zeta function ((■) and in case p = 0 we take nf=o^(- ■ •) = 1- 

Expansion (1.13) is the direct analogue of (1-12) for the Airy kernel determinant and completes part (A) 
of the aformentioned procedure for the same object. 


Remark 1.5. Note that in the limit ^ ^ -l-oo, all factors in the product (1.13) contribute to the 

error term and we correctly restore the well-known expansion for D{Jto leading order, cf. [ 1 , 19]. See 
Appendix A for further details on this matter. 


Remark 1 . 6 . Expansion (1.13) completely captures the behavior of D[Jai,i) in the underlying Stokes region, 
the Stokes curves in the {v,t)-plane are given by 

2 r 1 

■y = 3v2t-xlnt, X = 9-2’ 9eZ>i. 

In previous work [7] the very first Stokes curve q = 1 was identified and here we extend the results (with 
different techniques) to the full Stokes region. 

The most important consequence for us is contained in the following Corollary which is part (B) of the 
procedure for the Airy kernel. 


Corollary 1.7. For any fixed i G Z>o, we have, as s ^ — 00 , 

1 - A*(JAi) = ^25*+tf+3e-i^‘(l -t 0 ( 1 )), t = (-s)5. 
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This expansion matches exactly the formal result of Tracy and Widom in [36], (1.23). Next we turn our 
attention towards (1.4). 

Theorem 1.8. Given y S K and a > —1, determine p = p{x) ^ ^>o such that p = 0 for y < —4 and 
y+| <P<y+| for y > — 5 . There exist positive constants to = to(X) o) and vg = vo(x, a) such that 

D{JBess','y) = exp -^+01/5 (1 + ^1 + O ^max ^^ 

(1.14) 

uniformly for t = ^ tg^v = — ln(l — "f) > Vg and v >2t — 2(y + |) Int. Here 


di{a) = i\T{l + a + i)^ ^2 


— In—42—2a—3 


G(l + a) 

(271-) 5 


involving Barnes G-function G(-) and r(-) is the Euler gamma function. We take JliLo (■ ■ •) = 1 for p = 0. 


Remark 1.9. As j ^ + 00 , we rediscover the leading terms 0 /Zl( Jeess! 1) since all factors in the 

product of (1.14) move to the error term, compare [20, 25] as well as Appendix B. On the other hand the 
Stokes curves of D{JBess',x) ore given by 

p = 2t-2 (^y + Int, y = q.-^, g G Z>i 

and these have not been analyzed previously. 

The Corollary below summarizes the anticipated eigenvalue asymptotics for We have again match¬ 

ing with the formal result of [37], (1.27). 


Corollary 1.10. For any fixed i G Z>o and a > —1, we have, as s ^ -foo, 

^ 2ii+2a+3 


1 '^l('^Bess) — .. 


il r(l -I- a -f i) 


2*+l+ae-2t(i_^^(l)), i = 


Remark 1.11. With [36, 37], Theorem I .4 and 1.8 can be used to derive transition asymptotics for the 
underlying Painleve II and III transcendents. 


1.2. The ring of integrable integral operators. The technical part of the procedure is contained in its 
part (A), i.e. the derivation of transition asymptotics of, say, type (1.13) and (1.14). These asymptotics are 
obtained through an application of the Deift-Zhou nonlinear steepest descent method [21] to the following 
master Riemann-Hilbert problem (RHP). This RHP is associated with integral operators of type (1.10) and 
first appeared in [28]. 


Riemann-Hilbert Problem 1.12 (Master RHP). Determine Y{z) = Y{z;J,x) G a matrix-valued 

piecewise analytic function which is uniquely characterized by the following four properties. 

(1) Y = Y(z) is analytic for z G C\J and J C C is assumed to be a simple oriented curve. 

(2) The limiting values Y±{z) from either side of the contour J are square integrable and related via the 
jump condition 


f 1 — 2TTij(f(z)ip(z) 27rij<f^(z) \ 

— 27ri7'!/)^(z) 1-I-27ri7())(z)'0(z)y ’ 


z G J. 


(3) At possible finite endpoints of J, the function Y(z) is assumed to be square integrable. 

(4) As z ^ 00 , in a full neighborhood of infinity. 


Y{z) = I + Y^z-^ + O (z-2) , Too = 


For a given kernel (1.10) the derivation of an asymptotic solution for RHP 1.12 is kernel specific, never¬ 
theless the general philosophy is always to obtain first a “local identity” for the logarithmic derivative 


cHnZl(J; 7 ), (1-15) 

taken with respect to endpoints of J and 7 fixed. This means an identity for 9 In Z1(J; 7 ) involving local 
characteristica of F(z) such as its residue at z = 00 or behavior near specific points of J, see (2.8) and (6.4) 
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below. Second, through an application of the nonlinear steepest descent method [21], one derives asymptotic 
expansions for the local characteristica and, third, integrates (1-15) to obtain an expansion for D{J;^). 


Remark 1.13. The RHP 1.12 has been asymptotically analyzed for a variety of kernels over the past 20 
years. In particular in the context of gap probabilities D{J; 1) we mention, for instance, [17, 29, 15, 8]. Away 
from 7 = 1 the method has been successfully used in [9, 6] . 


1.3. Non-generic kernels in Hermitian matrix models. The Airy kernel (1.3) is obtained from scaling 
matrices in the GUE at the edge of the support of the eigenvalue densities p. Near a, say, right edge point 
X* of an interval in the density support we have 

p{x)c{x* — , xfx*, c > 0, fc S Z>o (1-16) 


and (1.3) simply corresponds to the generic situation fc = 0. In more fine tuned cases, i.e. for k € Z>i, the 
limiting kernels are described by a Lax pair solution associated with a distinguished solution of the (2fc)-th 
member of the Painleve I hierarchy. We will briefly discuss the first non-trivial case k = 1 and refer to [16] 
for a rigorous derivation of the underlying critical kernel: 


Kp2(A,p,;x,T) 


i ^n(A;x,T)^2i(p;x,T) - ^ii(p;x,t)^2i(A;x,t) 
2tt X — p, 


; Jpf 


(s, oo), s S K 


(1.17) 


with x,T gM. and the matrix entries ^jkiC) (viewed as analytic extensions from argf G (0, ^) to the entire 
complex plane) are characterized in terms of the following RHP. 


Riemann-Hilbert Problem 1.14. Determine $(C) = $(C; 2 :,r) G with x,r S M such that 

(1) $(C) is analytic for f G C\ 

ri = [ 0 ,oo). Pa = (-00,0], r2 = e“‘?^(-oo, 0 ], r4 = e‘T(-oo, 0 ] 

and the rays are oriented as shown in Figure 1.^ 

(2) The limiting values 4*±(C) from either side of the jump contour satisfy the jump relations 

$+(C) = <i>-(C)(-°iJ), CePa; $+(C) = ‘i>-(C)((, 0- C e Ti 

and 

d>+(C) = $-(C)([?), CGr2ur4. 

(3) $(C) is bounded at C, = 0. 

(4) As C —>■ oo, away from the jump contours 


$(C) = C-J^ 


p_(l 1 

C21-1 1 


I — ufTaC ^ + 


1 


2 -m\ 




c +0 


(ri) 


,^-6»(C;K,T)cr3 


where u,v do not depend on f, we define C,°‘ : C\(—oo,0] -G C such that > 0 for C > 0 and 

6'(C;a:,T) = - 2 xCT (1.18) 


Remark 1.15. The particular scaling chosen in (1.18) differs slightly from [16] and also [13]. For instance 
iI'(C;s,ti) in [13], Section 2, connects to $(^;x,t) in RHP l.lj via 

4 

$(C;x,t) = 2“tJ 5( , ( e c'y |J P^, x,tGR. 

1=1 


On the other hand (1.18) matches exactly [15], (1.21). 


The above RHP characterizes a solution of the second member of the Painleve I hierarchy, the Pf equation: 
indeed it is straightforward to show that $(^;x,r) solves a Lax system 


dx 


(C;x,t) = A(C;x,r)$(C;x,r), 


54) 


(C;x,r) = R(C;x,r)$(C;x,r) 


^Except for the opening angles of T 2 U r 4 the jump contours in RHP 1.14 and RHP 2.1 below are identical. 
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and its compatibility yields for u = u{x] t) and v = v{x] t), 

/5 5 1 \ 

(Pf) : X = -Tu - f -u^ “ ^ ^ (1-19) 

The precise expressions for and B{Q can be found in, say, [13] but they will not be relevant for us, we 
only use the existence of a real-valued, pole-free solution of (1.19) for x,t e M, see again [13]. Equivalently 
we only use solvability of RHP 1.14 for a;, r S K. 

Remark 1.16. The solution to RHP 1.14 is unique within the family of solutions of the form (j), )<i>(C; a:, r) 
with Lo independent off. However, the critical kernel Kp 2 {X, p,]x,T) is invariant under this gauge. 

Using the Lax equation ^ as well as symmetry constraints of RHP 1.14 we have for a;, r € K, 

i 1 f°° _ 

Kp 2 {X,p,;x,T) = — j ^ii{X;x + t,T)^ii{p;x + t,T)dt= — j $ii(A; a: -f t, x-f t, r) df, 

which should be viewed as the generalization of (5.9). The last identity implies in particular that Kp 2 : 
L^((s, oo); dA) O is positive-definite, compare Remark 1.3. For the analysis of the spectrum we would have 
to solve RHP 1.12 tailored to the choice 

4>{z) = $ 11 ( 2 ; X, r), tp{z) = X, r), Jp 2 = (s, 00 ) 

as s —)■ — 00,7 t 1 (in an appropriate Stokes regime) and x, r € M are kept fixed. The derivation of this 
solution is in fact very similar^ to the analysis of Kpj, presented in Sections 2 and 3 and it would provide 
us with an expansion for D(Jp 2 -^) of type (1.13). From it we would then obtain the desired information 
on the spectrum {Ai(Jp 2 )}i>o. We shall devote a separate publication to the spectrum analysis associated 
with non-generic edge behavior (1.16) for general k G Z>i. 

Remark 1.17. The purpose of this small section is to emphasize that the method of commuting differential 
operators will become very involved, if not impossible, when applied to, say, (1.17). As can be seen from 
the Lax system a commuting differential operator would necessarily involve a Painleve transcendent. And 
besides (1.17) several other kernels in the theory of random determinantal point processes (for instance 
Pearcey [38, 5] or Tacnode [22] kernels) fit naturally into the Riemann-Hilbert based scheme rather than the 
method of commuting differential operators. 

1.4. Outline of paper. The manuscript is split into the following two major parts. 

(i) First, in Sections 2 and 3, we derive an asymptotic solution to RHP (1.12) subject to (2.1) and (2.2). 

As can be seen in particular from Section 2 the steps we carry out in that section are closely related 
to the ones chosen in the analysis of the gap probability D{Jau 1), see [15], Section 3. The effect 
of 7 1, in contrast to the gap probability, becomes fully visible only in Section 3: we have a new 

g-function and we require different model functions than the ones chosen in [15], some differ only by 
a rank one perturbation, others are entirely new. In particular the new ones encode the discretized 
effect of the Stokes region and we use classical Hermite polynomials for our construction. 

Once all local contributions are in place we derive our first small norm estimates, however addi¬ 
tional steps are required to resolve a certain singular structure, see Subsection 3.5 for all details. The 
necessity of solving RHP 1.12 is related to the existence of differential identities for lniA( Jai; 7 ), com¬ 
pare Subsection 2.2. Using these identities we first obtain an asymptotic expansion for ^ ln_D( Jai; 7 ) 
and then perform a definite integration using the known expansion for Il( Jai; 1). In fact the known 
expansion for £>( Jai; 1 ) has to be improved to fit our purposes and we carry out the necessary steps 
in Appendix A. Finally, the information on {Ai(JAi)} is derived in Section 5 using an inductive 
argument. 

(ii) Second, in Sections 6 and 7, we focus on RHP 1.12 with (6.1) and (6.2) in the background. To the 
author’s knowledge the Bessel determinant D(JBess]"/) has not been analyzed previously with the 
help of RHP 1.12, although the analysis displays a few overall similarities with part (i). In fact we use 

"^This can already be seen from the analysis of the gap probability D(Jp 2 ', 1) presented in [15]. 
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again rank perturbations of more standard model functions and also here orthogonal polynomials, 
this time Laguerre polynomials. The asymptotic expansion for In Z?( Jeess! 7 ) is derived in Section 8 
and subsequently integrated using a refinement for £>( Jeess; 1) derived in Appendix B. The eigenvalue 
expansions are obtained via the same argument as in part (i), compare Section 9. 

Remark 1.18. In recent years several model functions in nonlinear steepest descent analysis with discretized 
parameter values have appeared, we mention in chronological order [14, 3, 4, 11]. These works have all in 
common that orthogonal polynomials (with the corresponding degrees as the discrete parameter) were used in 
the relevant constructions. 


2. Nonlinear steepest descent analysis associated with Kai - part 1 

The Airy kernel (1.3) is of type (1.10) with 

(({z) = Ai{z), -ifiz) = Ai'(z); JAi = (s, 00 ), sSK ( 2 . 1 ) 

and we shall derive an asymptotic solution of RHP 1.12 for sufficiently large (negative) s and 7 close to 1 
such that 

>fAi = ^ = T =-ln(l- 7 ) > 0 , t=(-s)i; X e M>o. ( 2 . 2 ) 


2.1. Preliminary transformations. Before we display the connection between Y(z) and D{J ai, 7 ) we first 
simplify RHP 1.12 in case of (2.1): Introduce the entire, unimodular function 


4>o(C) = y^e-'5"^(J 0;) 

and assemble an Airy-type parametrix, 

^(C) = ^o(C) < 


^Ai(C) e*3Ai(e 
^Ai'(C) e-‘t Ai'(e-'^ C), 


c e C 


(2.3) 


5 

argCe( 0 ,f) 


argC e (f^, 7 r) 


argCe(-f, 0 ) 


argC e (- 7 r,-f^) 


(2.4) 


This matrix valued function solves the following standard model problem which, in one form or another. 



Figure 1. The oriented jump contours for the Airy parametrix $(()) in the complex (^-plane. 
has appeared numerous times in nonlinear steepest descent literature. 

Riemann-Hilbert Problem 2.1. The Airy parametrix 4>((C) satisfies the properties below 
(1) 4>(C) is analytic for ( G T^- with 

ri = [ 0 , 00 ), Pa = (- 00 , 0 ], r2 = e“‘5(-oo,0], r4 = e‘?'(-oo, 0 ] 
and all rays are oriented “from left to right”, compare Figure 1. 
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(2) The limiting values from either side of the jump contours are related via 

$+(C) = $_(C)(_°iJ), Cera; $+(C) = $-(C)(Jl), CeTi; 
$+(C) = $-(0(1?), Cer2ur4. 

(3) $(C) is bounded at = 0. 

(4) C valid in a full neighborhood of infinity off the jump contours, 

1/1 1 

I e < 1 + 

48C 

where C“ is defined and analytic for f € C\(—c»,0] such that C“ > 0 for C > 0. 


= lie— 


1 /1 6i 

^ Ui -1 


O 


(rt) 




The model function (2.4) is useful as it allows us to reduce RHP 1.12 to a problem with constant, that is 
z-independent, jumps. This step appeared in [15] and amounts to the following transformation; We set for 
s < 0, see Figure 2, 

'T 


x{z) = y(z)$(z) 

and in case s > 0, see Figure 3, 

x(z) = y(z)$(z) 





z € u ^^2 u ^23 

Z € ^^4 

z G ^5 

z G U ^^2 U ^23 
z G ^4 
z G r25. 


(2.5) 


( 2 . 6 ) 



Figure 2. “Undressing” 
of RHP 1.12 in case s < 0. 
Jump contours of X{z) as 
solid lines. 



Figure 3. “Undressing” 
of RHP 1.12 in case s > 0. 
Jump contours of X{z) as 
solid lines. 


Hence we obtain for X(z) the RHP below. 

Riemann-Hilbert Problem 2.2. Determine X{z) G such that 
(1) X{z) is analytic for z G C\ Ui=i with 


t ( s ) 


= (s,oo), Pg®^ = (— 00 , s), = s + e ‘3 (— 00 , 0 ), F 4 *^ = s + e ‘3 (—c», 0) 


(2) The following jump conditions relate the limiting values X±(z), the jump contours are shown in 
Figures 2 and 3 as solid black lines: 

X+{z)=X.{z){^,l), zGVi^^-, X+{z) = X.{z){l°), zGTi^^uri^^ 

and 

x+(z) = x_(z)(i V)| 
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(3) In a full vicinity of z = s, 


Xiz)=Xiz){^^ 



arg(z - s) e (0, ’ 

1 0\ 

-11 j> 

arg(z - s) e (^ 

1 -1 W1 01 

0 1 All/’ 

arg(z - s) G (tt, 

1 -11 

0 1 /’ 

arg(z-s) G (^ 


27r > 
3 > 


Att ^ 


(2.7) 


with X(z) analytic at z = s and we fix the branch of the logarithm with arg(z — s) € (0, 27r). 
(4) As z ^ oo, 


X(z)=z-^‘ 


V2 V-1 


1 1 


1 ) { 


I + XcoZ ^ + XooZ + O 


(-*)} 




where we choose principal branches for all fractional exponents. The matrices XocX^o are z- 
independent, 

1 


oo 2 


i 1 




12 , 


X =11 ^ i 

°° 2 l-i 0 


(yii - r22) 

^■*^00 oo / 


2.2. Differential identity. The connection of D{JAi;i) to the solution of RHP 2.2 has been established 
for 7 = 1 in [15], (2.17). For 7^1 the derivation in loc. cit. can be copied almost verbatim and we simply 
state the following result. 


Proposition 2.3. For fixed 7 < 1, we have 

d 


-\nD{JAl.^) = ^^{X-Hz)X'{z)) 


21 


0 = - 
^ ’ Az 


( 2 . 8 ) 


in terms of the solution X(z) of RHP 2.2 and the limit is carried out for a.ig{z — s) S (0, ^) 


Remark 2.4. The requirement 7 < 1 is imposed for technical purposes only, it is in this case that RHP 2.2 
is solvable for sufficiently large v and t. For 7 > 1 the problem is only solvable for large {v,t) away from a 
discrete set in the {v,t)-plane. 

Remark 2.5. Another differential identity can be derived for lnZ)(JAi; 7 ), see Section A.3 for further 
details. 


On the upcoming pages we will derive an asymptotic solution of RHP 2.2 and then compute the asymp¬ 
totics of ln£)(JAi; 7 ) through ( 2 . 8 ). 


3. Nonlinear steepest descent analysis associated with Kai - part 2 
3.1. Initial transformation. From now on we assume that s < 0 is sufficiently large negative. Define 

T{z)=X{\s\z + s), zeC\(SrU{0}). 

which “centers” the problem at the origin z = 0 so that we have jumps on the contour 

4 

St = U r, 

1=1 

shown in Figure 1. More precisely we obtain 

Riemann-Hilbert Problem 3.1. Determine a function T{z) = T{z;s,^) G which is uniquely char¬ 

acterized by the following properties: 

( 1 ) T{z) is analytic for z G C\(St U { 0 }) 

(2) We have the jump conditions 

T+(z) = r_(z)(i V)- ^eFAW; T+(z) = r_(z)(AJ), ^eraUO} 

and 

r+(z) = r_(z)(]o), ze (F^uFAUO}. 
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(3) In a neighborhood of z = 0, 


T{z)=T{z){l 


arg0e(O, 

ii?), argze(^,7r) 

(Jli)(} 0 ), argze(^,f) 
argz e (f^,27r) 


where T{z) is analytic at z = 0 and the branch of the logarithm is specified by the requirement 
arg2; € (0, 27r). 

(4) The normalization at z = oo reads as 


Tiz) = {\s\z) })e-'f‘^={/ + Xo.(|s|z)-5+0(z-i)} 


e-|(k|z+s)2o'3_ 


Remark 3.2. Compared to [15], the RHPs 2.2, resp. 3.1 also have a jump on (s,oo), resp. (0, c») since 
we are interested in the analysis of with 7 ^ 1 m general. This will have a crucial impact on the 

steps below. 


3.2. Normalization transformation. We define for z S 


g{z) = -z^ [z- -] +Vln 


1 + (2z)2 

1 - (2z)5 




(3.1) 


with principal branches for all fractional exponents and logarithms. In particular we choose 

l-(2z)i) - 

Further steps require the following analytical properties of the 5 -function: 


—TT < arg 


Proposition 3.3. The function g{z) introduced in (3.1) is analytic for z G C\((— 00 ,0) U (|,-l-oo)). In more 
detail along the real axis with orientation as shown in Figure 1, 

9±(--) s ± ie) = ±?iv7I (" - |) + ‘'""8 ( ]^]([|| ) ' ^ 


well 


and 


Also, 


9±(z) = (^z- 0 -f Pin 


1 - 
l' 


, z G 0 


y2z- 1 


± iirP, z G ( -, -boo 


near z = oo, 


—-(z - 1 ) 2 - 1 - g{z) = iirP + ^V\/2 — z ^ + O [z ^) , z ^ oo, z ^ 
At this point we introduce 

S{z) = z G C\(St U {0}), 

which leads us to the problem below. 


Riemann-Hilbert Problem 3.4. The normalized function S{z) = S{z;s,'y) G is characterized by 

the following properties 

(1) S{z) is analytic for z G C\(St U {0}). 
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(2) The limiting values S±[z),z G Et are related by the equations 


S+{z) 

= S-{z) 1 

^ 1 0 \ 
^e 2 ts( 2 ) 1 J 

and 




5'+(z) 

= S_{z) 


,S+(z) 

= S-{z) 


2 g-t(>CAi+ 2 s( 2 )) 

0 1 


’ ^ ( 0 ; 2 ' ’ 


i‘( 9 +( 2 )- 9 -(z)) Q-t{^Ai+g+{z)+g-{z)) 


0 


-t(g+( 2 )-s_( 2 )) )> ^elx,+00)- 


(3) Near z = 0, with argz G (0,27r), 


ei7rtVa3 5.(^)e-t9(2=)a3 ^ J ^ 2,ri^ 


(4) have the normalized behavior at z = oo ^ Et, 



argzG ( 0 ,^) 

1 01 
-11 j: 

argz G (f^, 7 r) 

1 -nn 01 

0 1 All/’ 

argzG (tt, ^) 

1 -i\ 

0 1 /’ 

argz G (^,27r) 


5(z) = |/+ + o (^-i)|, 


with 




We make a few important observations: With 0 < r < | fixed, 

3^(5(-2)) < 0, z G (r 2 U r 4 )\i:)( 0 ,r); L>(zo,r) = {z G C : |z-zo|<r}; 


(3.2) 


next for z G (0, |)\(I?(0, r) UD(i, r)) and sufficiently large t > to,v > vq (see (2.2) for the definition of 


., + 2sL) = X., + ivL - 0 + 2V 1„ j i 


> 0 . 


Finally, 


g+{z) — g-{z) = 27riF = const. G iM, G ( 2 > ) ’ 


(3.3) 


(3.4) 


+ 9+iz) + 9-iz) = + +‘2.V\n{^ 


y/Tz + l 

v^- 1 


> <5 > 0, 2 G -,+oo \D -,r . 


1 


These estimates lead us to the expectation that the major contribution to the asymptotic solution of the 
S'-RHP arises from the line segments (—oo, 0) U (|, +oo) as well as two small vicinities of z = 0 and ^ = a. 


Remark 3.5. Observe that through (3.1), 
and if we represent 


tv = X, 


thus 


X = k + a; A: G Z>o, 


^TTitV _ g27ria 


We now continue with the relevant local analysis. 
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3.3. Analysis of model Riemann-Hilbert problems. The outer model function^ 

z € c\(^{-oo,0] U [^.+oo)) 

with the scalar Szego function 

'i + ( 2 z)A" 


(3.5) 


Viz) = 


z G 


2 a 


P(z)=e-“(l + ^+0(z-)), 


1 - ( 2 z )2 J 
satisfies the properties listed below. 

Riemann-Hilbert Problem 3.6. The parametrix P^°°\z) has the following analytical properties 

(1) P(“)(z) is analytic for z G C\((— oo,0] U [i,+oo)) with orientation of the real axis as indicated in 
Figure 1. 

(2) The function P^°°l (z) assumes square integrable limiting values on (— oo,0] U [^,+oo) which are 
related by the jump conditions 

p|“)(z) = J), zG(-cx), 0 ) 


J_), zg{^,+oo 


(3) As z —>■ oo, z ^ K, 

P(“)(z) = (|5|z)-J‘^3e-™^ e-' 


4<73^17ratT3 


{7 + o(,.-i)}. 


The local parametrix near z = 0 differs from the one used in [15], Section 3.4, by a rank one perturbation, 
compare (3.10) below. In more detail we first define 


-ciKMhJ- 


CgC\(-oo,0] 


(3.6) 


in terms of the modified Bessel functions Iq and Kq and with principal branches for : argC G (— 7 r, 7 r]. 
The standard properties of Bessel functions [33] in mind we obtain a bare Bessel parametrix: 

Riemann-Hilbert Problem 3.7. The function J{f) G defined in (3.6) has the following properties 

(1) J(C) is analytic for C, G C\(—oo,0]. 

(2) On the negative half ray, oriented from —oo to the origin, we have 

A A 

,0 1 , 


J+(C) = J-(C) 


(3) Ts C ^ 0, 


^(C) = J(C)(J 


, Cg(-oo,0). 

, argC G (- 7 r, 7 r), 


(3.7) 

(3.8) 


and J(C) is analytic at f = 0. In more detail. 


J((j) = e A'^3^2'^3 


/o(C^) 


i(/o(C5)ln2 + Er=oV'(fc+l)®^) 


- (c^^o(C^) ln2 - loiC^) + J2T=i + 1) 


ICI< 


with 


oo (l^\k 

= d7;(d) = E“,t!)7 


(-O'" r'idi 

f:{C) = ^, CgC\{ 0 ,- 1 ,- 2 ,...}. 


k=0 ' ' k=l 

(4) The function J{C,) is normalized so that 


no 


1 




1 1 


1 /-I - 2 i 


8C3 


- 2 i 1 


+ 0{C^) 


as f ^ oo with —tt + (5 < arg C, At: — 5 and (5 > 0 fixed. 


(3.9) 
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We can now define the parametrix near the origin in terms of the model function J(C): For 

0 < |z| < i 


gi 7 rtVa 3 p( 0 )(^)g-ts(z)a 3 ^ E^^'> {z) J {(:{z)) ( J 


I, argz€(-^,^) 

(_\?), argz€(^,T) 
( 10 ), argz€(-^,-f) 


with the locally analytic left multiplier 

£;(o)(z) |z| < I, 

8 

the locally conformal change of coordinates, for \z\ < 


C(^) = t { giz) - - In 


a (1 + {2z)i 


t Ii-(2z)5 




1 - - 


41 + 2if 


3l-2if 




(3.10) 


and the branch of the logarithm in (3.10) such that argz G (0,27r). It is straightforward to verify the 
analytical properties of P*^o)(^ 2 ;): 


Riemann-Hilbert Problem 3.8. The parametrix P^o) ( 2 ;) has the following analytical properties 

(1) P^‘^\z) is analytic for z € D(0, |)\(I]t U {0}) with D{zo, r) = {z G C : \z — zo\ < r}. 

(2) Since f = f{z) locally conformal near z = 0, we obtain directly the jump behavior 


pf(z) = Pl°^(z)(g2i.) J), ze ((r2ur4)\{0})np(^0,i); 

pf(z) = p!°^(^)(^_°^ J), z e (FaUO}) n p (^0, 


compare Figure 1 for orientation, and where we used (3.7) in the last identity. Also, by the choice 
of branches in (3.9), 


P®(z) = Pjf’{z) 


1 ( 0 )/ 


^-t(>CAi+2g(z)) 

1 


ze (ri\{0})np 0 



All together, P^^\z) models precisely the jump behavior of S{z) for z G D(0, |)\{0}, see RHP S.j- 

(3) Near z = 0 with argz G (0, 27r), we deduce from (3.8) and (3.10), 

(I, argze(0, f^) 

iii-tVCT3p(0)/ \ -i3(2)o'3 _ p(0)/ I A 2^ J (-1 1 )’ argzG(^,7r) 

1 / argzG(.,f) 

UoA)> argz G (f^,27r) 

which matches exactly the singular behavior of S{z) near z = 0. 

(4) As t ^ + 00,7 t 1 subject to (2.2), we derive from (3.9), 

p(0)(^) ^ p(oo)(^) |j+ _^(p(^))-3 -2i^ (P(z))"^ +0 (t-2)| (3.11) 

uniformly for 0 < ri < |z| < r 2 < |. Here, we used in particular that on the latter annulus, 

^-t{^A,+gPA+a-(A) = o {t-'^) , t^+oo,7tl: x^, = lV2-x—. 

3 i 


Remark 3.9. Note that from properties (2) and (3) in RHP 3.8 we obtain 

Siz) = No{z)P^°\z), 0<|z|<i 


where Nq{z) is analytic at z = 0. 
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The outstanding parametrix near 2 ; = ^ is in some sense more elementary than the Bessel-type parametrix 
(3.10). We draw inspiration from [14, 3, 4, 11] and define first 

PkiO 


H{0 = 


Tt-, f rii ® ’ CgC\k, fcez>o 

Ik-iPk-iiC) t=c/ 


(3.12) 


27ri JkPK-iUI'" t-C 

with the help of monic Hermite polynomials {pfc(C)}feGZ>o:P-i(C) = 0 = 7 _i, cf. [33]: 

Pfc(C) = C" + afe.fc-22'=-' + 0(C"-"), C 

where 


/ Pj{t)Pk{t)e ^^dt = hkSjk, 

Jr 


x/t 1 

hk = k\^, fc G Z>o; ak,k -2 =-^k(k - 1), kGZ>o. 

Remark 3.10. The normalization in [33] of the Hermite polynomials {R/c(C)}feGZ>o different from the 
one chosen here: we have to use the relation pk{C) = H}.{C,),C, G C. 

These properties lead at once to a bare Hermite parametrix: 

Riemann-Hilbert Problem 3.11. For any k G Z>o, the function H{f) G defined in (3.12) has the 

following properties 

(1) i?(C) is analytic for C, G C\K and we orient the real axis from —00 to -l-oo. 

(2) Along the real line we have 

A A 

.0 1 , 


R+(C) = R-(C) 


CG 


(3) Ts C “t 00 ,^ ^ R, with 7 fc = — 


1 - 


Ik 

0 


0 

7fc-l 
7fe + l , 


0 lkl2\ 

(3 V“7fc-3 0 ) 


+0 (C^) 


For the actual model function we then take with \z — 


where 


p(5)(z) = ^;(5)(2)i^(((z))e5^^‘'"3e‘s(^)‘"^ 


(3.13) 

(3.14) 


< 


"-2 


< 8 ’ 


with the choice 


m = az) 


l + {2z)i 
1 - (2z)5 


X , ^ \ h ^ 5fc o: I ^ 2 /l / j. 

t 2 =(-l )'=2 4t 2 Ji+ Iz- 


O 


^ 2 


is analytic at z = \. Throughout, ^ = C( 2 ),l 2 :—|1<| denotes the locally conformal change of variables 




and we use again the representation 




^>0 9x = fc-l-Q!, k & Z>o, —- < o < -. 


2 - 2 

The important properties of pA)(^z) are summarized below. 

Riemann-Hilbert Problem 3.12. The parametrix pAl (^z) has the following analytical properties 
(1) pA)(^z) is analytic for z G D{^, |)\Er. 


(3.15) 
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(2) By local analyticity of ( = C('Z)) we have 

'X 0-t(>^Ai+2s(2;)) 

.0 1 


= P'x}\z) 




= P^}\z) 


g27ritV g-t(>fAi+S+(z)+S-(2)) 


0 


^—27TitV 


( ^,+00 ) ■ 


This matches exactly the jump behavior of S{z) for z G compare RHP S.j. 


(3) As t ^ + 00,7 t 1 subject to (2.2), we derive from (3.13), 


p(i)(^) = p(-)(^)<)/+_(p(^))- 


7fc_i/32(z) 


(!>(=))" + 


+ 


0 (z) 




lk+2 P (^) 




^-7fc_3^2(z) 0 

uniformly for 0 < ri < \z — ^\ < r 2 < ^ ■ 

Remark 3.13. Note that from property (2) in RHP 3.12 we obtain 


O 


CHz) 


7fc-2 

Ik 

0 


0 

7fc-i 
7fc+i ^ 


(3.16) 


S{z) = Ni{z)P<^^'>{z), 0< Z-- < 

where Ni{z) is analytic at z = 

This completes the construction of local model functions, we now use the explicit functions P^°°') (z), (z) 

and P'^^'>{z) and compare them to the unknown S{z) in RHP 3.4. 

3.4. Ratio transformation and first small norm estimate. With (3.5), (3.10) and (3.14) this steps 
amounts to the transformation 

, -1 


/I n\ 

As{z)l{p(h){z))-\ |z-i|<r ; u; = -\s\i{N{S^-aV2a:,)N-^) 

^ ^ [(p(-)(z))-\ |z|>r,|z-i|>r 

in which 0 < r < | is kept fixed. Here we have made use of the abbreviation 


21 


N = e 


1 


1' 


g-i^ (73gi7racr3 


v7V-l ly 

and Sao occurred in RHP 3.4. Recalling RHP 3.6, 3.8 and 3.12 we are lead to the following problem. 



Figure 4. The oriented jump contours for the ratio function R{z) in the complex z-plane. 

Riemann-Hilbert Problem 3.14. Determine R{z) = R{z;s,'y) G such that 

(1) R{z) is analytic for z € C\Sji with square integrable boundary values on the contour 

E^ = dP(0,r)uapQ, 

which is shown in Figure 4- 
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(2) On the contour we have jumps R+{z) = R-{z)Gr{z] 5 , 7 ), 2 : £ with 
GR{z;s,^) = P<^^\z){P^°°Hz)y\ zGdD{0,r); Gr{z- s,^) = {z){P<^^'>{z))-\ zea 7 ^Q,r^; 

Gr{z-,s,^) = P^°°\z)(^^ ^ J{P^°°\z)) \ 


{r,2~^ I ’ 


GR{z;s,^) = Py\z) 


1 g2'n-iae-t(>CAi+g+(2)+g-(i)) 

0 1 


yy\z)) , ze[-+r,+oo 


Gr{z;s,j) = P^^\z)(^J^,) ^^{p^°°Hz)) \ ze (r2ur4)n{zeC: |2| >r}. 

By construction, see Remarks 3.9 and 3.13, there are no jumps inside D{0,r) U Dy,r) and on 
{— 00 ,—r). Moreover R{z) is bounded at z = 0. 

(3) As z ^ 00 , 

^ (<i 1) {/ + 5 'oo 2-5 +c>(2-i)| 

= / + 0 ( 2 - 5 ) . 

Through standard small norm estimations, compare (3.2), (3.3), (3.4) and (3.11), we obtain at once 

Proposition 3.15. Given x € I^>o there exist positive to = to{x),vo = vo{x) c = c{x) such that 

c 2 

||Gfl(-;s, 7 ) - /||L5nL<»(E„\aD(i,r)) < jr, yt>to,v = -ln(l - 7 ) > vq ■■ v= - - xlnt. 

The circle boundary dDy,r) requires further analysis: From (3.16), uniformly for 2 £ dD{^,r), 


Gr{z;s,x)-I = (|s|2)-i'^^{p+(2)t-H« + p-(^)t-|-«+P2(^)i-i+p+(^)t 

+ p^(2)ri-“ + 0 (t- 2 +i“i) }(|s|2)^‘"^ 


“ 9+“ 


(3.17) 


with the t-independent matrices 

p+(z) = -i7-i^J^e-*-' 

^ 2 C( 2 ) 

which are both of rank one, and 


1 1 
-1 -1 


R^iz) = 

2({z) 


1 -1 

1 -1 


R2iz) = 


ACiz) 


as well as 


1 k 
k 1 


R 3 (^) = -m-3 


R3 (z) = -Gk+2 




2C3(2) 


1 1 

-1 -1 




iTratTs 


1 -1 
1 -1 


iiracTa 


2(3(2) 

We have introduced j8{z) = t~^^z) and (( 2 ) = t^^z) where / 3 ( 2 ) and (( 2 ) are t-independent. Note that 
for a £ [— 5 , 5 ) the first two leading terms in (3.17) are in general not close to zero. In order to overcome 


(3.18) 
1 


this feature we use matrix factorizations, 

Gr{z-s,x) = E+{z){\s\z)--^^-[l + R^{z)t-^^-^ + {R^{z)-Rt{z)R-yz))t-^ 

+ (P+(2) - P+(2)P2(^))t-®+“ + R3{z)t-i~" + O (t-2+l“l) }(|s|2)i'^3^ 0 < a < i; 
Gr{z-,s,x) = E-{z){\s\z)-z-^[l + R+{z)t-"^+^ + {R2{z)-R^{z)Rt{z))t-^ (3.19) 

+ (P3-(2) - R-yz)R^{z))t-i-^ + P3+(2)t-i+“ + 0 (r2+l“l) }(|s|2)i-^ -^ < a < 0, 
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where 


are invertible and meromorphic for z G D(^,r). The factorizations (3.18), (3.19) above, i.e. 

Gji(z;s,j) = E^(z)E^(z), zGdnQ,r 

motivate our next move. 


3.5. Singular Riemann-Hilbert problem and iterative solution. In this step we introduce 

ffE+(z), |z-i|<r 


Q(z) = R(z) < 


|z- 2 I > r 


0 < a < i 




|z- 5 I > r 


which leads us to a singular RHP. 

Riemann-Hilbert Problem 3.16. Determine Q{z) = Q{z;s,'y) G such that 

(1) Q{z) is analytic for z G C\(I]h U {|}) with square integrable boundary values on the jump contour 
T,ji shown in Figure j. 

(2) The jumps on are as follows, 

Q+{z) = Q-{z)GR{z;s,-f), z G Si?,\dZl Q,r^ ; 


Q+(z) = Q-{z) 


(E+{z), 0 < a < i 

\eo{z), -5<a<0’ 


z GdD [ -,r\ . 


(3) The function Q{z) has a first order pole at z = In more detail, near z = I, 


Q{z) = Q{z) 


vO 1 




^ 2 


< r; 


where Q{z) is analytic at z = ^ and we have introduced 


(T^ = 


. -lr,-SL_5 1 

— 17j^, 2 2 2 

1 ■ -lr,-fiL_S ,Q,_ 1 ’ 

1 — 17 ^, 2 2 4 2 


a G 




17fc l2 2 4t 


-I • r\ _H I _ 

l + i7fc_i2 2 it 


—, 01 G 


- 2 ’« 


as well as 
T+ = ( 


(4) As 


_e27ria(M)2 Q 


CX) 



l\ 

^ / 1 1\ 


1 

a G 

0 , - ; 


CJ G 

—,0 

’ 2j ’ 

le27na(l£i)2 q) 

2 


Q(z) =/ + 0 ( 2 - 5 ) . 

Remark 3.17. The singular strueture (3.20), (3.21), (3.22) follows from the observation that 

1 


R{z) = Qiz){E^{z)) \ 


"-2 


< r 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


is analytic at z = Hence the singular part of Q{z) can he derived by comparison in (3.23) which leads to 
(3.20). 
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Note that all jump matrices in the Q-RHP are close to unity at the cost of an isolated singularity at 
z = ^. This will now be resolved by a final transformation. Define L{z), z G such that 

Q(z)=-l^(^z-^^I + B^jL(z)^, (3.24) 

with G constant in z. 


Riemann-Hilbert Problem 3.18. Determine L{z) = L{z;s,j) G such that 

(1) L{z) is analytic for z G with square integrable boundary values on the jump contour shown 

in Figure 4- 

(2) The jumps are identical to the ones in the previous Q-RHP 3.16, i.e. 

L+(z) = L_(z)Gq(z;s,7), z G 

(3) The function L{z) is analytic at z = ^ provided 



which follows directly from (3.24) and (3.20). 

(4) As z ^ oo, we have that L{z) -G I. 



(3.25) 


Remark 3.19. Note that 







and thus 

tr 5^ = 0 = det . 

This implies that det((z — ^)I + B^) = (z — ^)^ and hence (3.24) is consistent with the identities 

detL(z) = det( 5 (z) = 1 . 


At this point it is clear that the L-RHP admits direct asymptotic analysis, indeed we have 

Proposition 3.20. Given x = k + a G K>o with k G Z>o, — 5 < a < 5 , there exist positive to = to(x): vq = 
Vo{x) and c = c(x) such that 

liG'Q(-;s,7) -^llL2nL“(aD(i,r)) < ct"s-l“l, Vt>to, u = - ln(l - 7 ) > iiq : v =‘^y/2t - x^^t. 


Hence together with Proposition 3.15, by standard arguments [21], the singular integral equation 

L{z)=I+^( L_{w){Gq{w)-I)^^, zGC\Er 

J'Sfi w — z 

which is equivalent to RHP 3.18 is solvable for sufficiently large t > tQ,v > vq. In fact 


(3.26) 


Proposition 3.21. Given x S ffi>o there exist positive to = to(x) 5 ^o = vo(x) c = c{x) such that RHP 
3.18 is uniquely solvable for t > to,v > vo ■ v = ^V^t — ylnt. The solution L = L{z; s, 7 ) satisfies 

[[-^-(ss^t) --f||L2(EH) < yt>to,v>vo. V =‘^V2t-x'^nt. 

The last Proposition can be used to derive an asymptotic expansion for the logarithmic s-derivative 
through (2.8). In the derivation of such an expansion a certain structural information will prove useful: 

Note that all jump matrices in the Q-RHP are in fact (formally) conjugated by compare for 

instance (3.18), (3.19), i.e. we can write 

Gq{z-,s,x) = |s^^'""GQ(z;s, 7 )|s| 3 '"^ z £ 'Br. 

This leads to a RHP for the function L{z) = L(z;s,j) = jsj i‘^3L(z; s, 7 )|s|“ which is uniquely solvable 
by iteration, 


li^^-(-;s,7) -< ct 2 l“l^ yt>to,v>vo: 


2 R- 

V = - v2t — ylnt. 
o 


(3.27) 
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We recall (2.8), 


4. Extraction of large gap asymptotics 


= ^lnZ?(JAi;7) = 

2tv ^ 


where the limit is carried out for arg(z — s) S (0, 

4.1. Asymptotics for the derivative. Through the sequence of transformations 

X{z) I—?> T{z) I—>■ S{z) i-A R{z) i-A Q{z) I—>■ L{z), 

identity (2.8) leads us to 

^ -{T-\w)T'iw)). 




27ri|s| 

7 

27ri|i 


>21 


ryQ 2i5^.(0) 

= — .. I — {S~ {w)S'{w)) 

W^o 27 ri|s| ^ I I V 7 


21 


Id—>-0 


{(P(0)(u;)) '(P(°^(ir))'}^^ 


7 


lu-s-o 27ri|s| 


{(P(°)(u.)) ^Q-1 (u;)Q'(r;)P(°)(r;)}^^ 


iti—>-0 


and the first summand is computed with the help of RHP 3.8 as 


7 


27ri|s| 

For the second, with (3.24), 




21 


= t7|s|' 
4 


1 - 22 - 
t 


7 


27ri|i 


27 


Id—>-0 


{(p(°) {w)y^Q-\w)Q'{w)P^°^ {w )} 

+ ^ {(p(°)(ic))-'T-'(ic)T'(ic)P(°n«^)} 

{l-1(0)P±P(0)} 


{(P(°) (w)) ^L-\w)B^L{w)P'-°^ (ic) 


It ;—>-0 


21 


IP—>^0 


= 270 


— 27riQ: 


1-22 


_ T —27ria 

21 2 


1-22- 


t 


[l-\q)l' 


We now begin to compute the contributions from L=*=^(0), T'(0) and L*^(i), T'(|). In order to achieve this, 
we use the integral equation corresponding to L{z). Modulo exponentially small contributions, for z G C\E/j, 

dw 


L{z)=I+^<f {dQ{w)-l)— + ^<f {dQ{w)-l)PPX_ (4.1) 

2711790(0,7-) 'w-z 27117a£,(i 'w-z 

+ ^(f (L.{w)-l){GQ{w)-l)— + ^(f {L_{w)-l){dQ{w)-l) — 

271-1 /aofo.r-l 'w-z 2TnJgiy(iy ' W - 


With (3.27), 


(L_(ic)-/)(Gq(r;)-J)—, zGC\Eo. 
2711790(0,7--) 'w-z V ) 


ldD(0,r) ' ' ' W-Z 

For an analogous estimation on dDy,r) which we require in (4.1), we use 


1 

27ri 


<f (P_(r;)-J)(Gq(r;)-/)—= o(t-i-2|“l), z G C\En; 

JdD(^,r) W-Z \ / 


and therefore back in (4.1), for z G C\So, since 1 < 1 + 2|a| < f + lal < 2, 

Hz) = I+7^f {GQiw)-l)^^ + ^ (f (GQ(ii;)-/)^^+G(t-i-2|“l). 

2711/90(0,7)^ w-z 2711 /a£,71,7 'w-z V J 


Starting from (4.2) it is now straightforward to derive the following estimates, 

|s|i<^3 7^-l(0)|5|-i<-3 ^7_e-i7ra<73|^±^-i-|a| |'^-min{l+2|a|.f-|a|}^ 


(4.2) 


(4.3) 

(4.4) 
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with 


and 


Aln = 1 - 2 !- 


-1 


0 


2 1 + 2 ! I 


Al, = -^ 1 - 2 ! 


,-3 0 

0 


1 + 25 


1-22 $ 


= 2-!-(fc± 1)2J'"3 
o 


Af = 2-!-(A:± 1)21'^=' 


0 ..3, 

1 o' 


0 13 
7 0 


2 “J" 


Remark 4.1. Explicit expressions for the matrices A^ and A 2 , for instance 

‘ ' . Tfe I (~) 


are not going to be important later on: As we shall see shortly, all powers of t which explicitly contain the 
parameter a will not contribute to leading orders after integration, see (4.10) and (4.12) below. 


We now obtain 


{l-\ 0)L'(0)}^^ = e2’^'“|s|!{A±t-!-l“l + A2t-^ + Aft-^ + O (^t-™{i+ 2 |«l.f-|«|}j 


' 21 

= 


(4.5) 


_1 + ± 1) ± 2±f+ir!-i“i J 

4t 6t ^ ’ 


i7fc^ 


(-) 




fl\ . , , 1 

/ 

1 1 1 \ ■ 

- = s - 1 '"= 

/T0 f 

t-!-|“l ) 

\2 ' ' 

V 

yj 


Since furthermore 

L±i 

we also deduce from (3.25), 

=CT±e-'’"“'"3|srJ'"3{/ + C>(r!-l“l) }2J 
and thus 

{l- 1 (o)b±l(o)} 


U|3^3eivra<T3 


1 ±1 
Tl - 


= e^™|s| 


21 


!(T±{23'^3 


1 ±1 

Tl -1 




2-4'"3 Q 


(4.6) 


Summarizing (recall that 7 = 1 — e = 1 T O (t °°)) 

2 


l^(^Ai;7) = 


- Ip (-1 + ± d) +0 (i-Ui‘i) 


4 8 s 

which is the central estimation for the upcoming integration. 


- ^ T V^|s|^cr^ - - — T :^(A; ± 1) T 0 (t ^ l“l'j 

85 s 12s V / 


(4.7) 


4.2. Integration of expansion (4.7). Let us first work out the details in the special case k = 0 and 


0 < a < Note that in this situation 


r(JAL;7) “ + “ ^ “ y2|s|++ + o 


(rS-), 


cr+ = 


7r-!2-tt“-! 

1 + 7r“5 2“4f““! ’ 


(4.8) 


and thus the following Lemma will be useful. 

Lemma 4.2. Let t = t{s) = (—s)! and sq < s < 0. For any f G L^{so, s), we have 
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Recall that 


so that 


>fAi =-^ = X —; T =-ln(l - 7 ) > 0 , x = k + a] 




(4.9) 

(4.10) 


and thus back in (4.8), 


—/ (—u)2a~^(t{u))du=-'/2 
J Sq 


' TT 22 4w 20 


V2i 


io 


1 + TT 22 iw 2e~'"e3'^' 


-dw 


= ln(^l+7r 22 iw 20 


^ 1 /■* TT 22 410 20 ^03'^'^ dw 


W=to 


/to 1 + TT 2 2 4UI 2 0-«0 3 


V 21 




Let us now choose {v,t) sufficiently large positive such that 


3 3 

V < t < 


v+-lnt] , 


2 V 2 ~ ~ 2 V 2 

i.e. the base point of integration is Sq = ~(^^)^ O'" equivalently to = Then 

—V 2 ( {—u)^a^(t{u))du = In fl + 7 r“ 52 “tt“ 5 ei'/^*“'“^ + O ■ 

J Sq 

Summarizing, we can integrate (4.8) as follows, 

/det(J-7i^Ai)L2(,^^) 


In 

\^det(/-7ifAi)|i2(,„,, 
where we used that 


= —(s^ — Sg)-In 

12 ^ 8 


+ ln(^l + 7r 22 it 20^3/2* 2 ^^ (4.11) 


f O (t-5-“(u)) du= f O dw= f O dw = O fri) , (4.12) 

«/ Sn ^ tn J to 


uniformly as t ^ + 00,7 t 1 such that v <t < | Int). However, with (A.l), we have 


that 


lndet(/ - 7 A:Ai)|^ 2 ( 5 ^^^) = + Incg + O (tg , 


Co = exp 


-In2 + C'(-1) 


since for the latter determinant > fv^- This substituted back into (4.11) we obtain 

i i 

= ^ - bn|s|+lncg + ln fl + T 2 -lr 563^*"^ +o(t-^) 

12 8 \ yTT J V / 


i(so) 3 

Proposition 4.3. There exist to > 0 and wq > 0 such that 

.3 

lndet(/ — 7 iLAi) 


uniformly for 


L2 {s,Oo) 


t > to^V > Vq ■ 


3 3 

v<t< 


2 V 2 


2 V 2 


u + 2 ^ 


(4.13) 


Remark 4.4. Note that the result of Proposition 4-3 matches the leading order expansion derived in [7], 
Corollary 1.17. We only have to use that (4.13) implies 


t>to, V >vo : 


2 /- V 2/- Int If fhrt 

-V2 >- = }<,,> -V2 - Y^ with Y = - 1 + O ^ 


t 


For the general case k S Z>g we require the following notation: Let 
3 


tfc — 


2 V 2 


Int^ , i'^. = + k\nt), k G Z>i; to =/'o 


= tg = u. 


or equivalently Sk = , s'f, = —iff. In this case (4.11) generalizes to the following two expansions 
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Lemma 4.5. For ik <t < i'j,, i.e. a S [—^,0], 

/ det(/— 7i^Ai)|^2(s Qo) \ 1/3 .. 3 \ 1, S , 


, / ^ _ 5^ _5\ Ik . - N , 

- In (1 + I 7 fc_i 2 2 4)+—(fc-l)ln 


s 

- 2fc2 In 

s 





+ —'/2k(t — ik) + O ^ 2^ _ (4-14) 


■) 

Proof. We use that 

V 2 f {—u)^a~(t{u))du = —-'/2 f 

Jsk 37 1 , 

= In ^1 + - In ^1 + i 7 /c_i 2 ^“J^ + O 

and the stated identity follows from (4.7). 


'* i7fe_i2¥-iii;'=-5e^e-i^“ 

ik 1 + i7fc_i2T“f e”?^’ 


■ dw 


□ 


Lemma 4.6. For i'^.<t< ik+i, *-e. a. e [0, ^], 

Uet(/- 7 AA,)L.„;„, ' 12' ” 


7k 

+ Y^(fc + l)ln 


S 

•'/ 

- 2fc2 In 

s 

'■f 

^k 


^k 


+ ln(^l-i7^^2 ^2*" 4i 2'^ 

+ 2'\/2fc(t — ffc) + O ^) ■ (4-15) 


Proof. Use 


-^2 / (-u)2a+(f(u))du=-72 ^ _fc_i ^ 2^^ du; 

Js^ o Jt'^ 1 — 17 ^ 2 2 4w ^ 2e“’'e3X2:i« 

= In — i7“7“7-ft-^-5e3^*“'"^ + C> > 


back in (4.7). 


□ 


The idea at this point is to successively improve Proposition 4.3 with the help of Lemma 4.5 and 4.6. For 
instance, start with fc = 1 in Lemma 4.5, then with Proposition 4.3, 


lndet(/ — yi^Ai) 
and back into (4.14), 
lndet(/ — yT^Ai) 


(si ,c 


12 8 


- In |si| + Inco + In + 71 ^2 


L^{s,oc) 12 8 


- In |s| + Inco + In ^1 + 7722^20 3 ^*+’'^ — In ^ 7122 ^^ 


=-In I s 

12 8 ' 


+-V 2 (t — ii) + O (^t 2 ^ 
where we used that 

^72(t-U) =ln (r^et^*-") . 

Summarizing, 

Proposition 4.7. There exists tg > 0 and uq > 0 such that 


l + TT 22 4t 2 ^ 


lndet(/ — yi^Ai) 
uniformly for 


1 


= - X ln|s| + Inco + In ( 1 + ^ 2 U ’ 

L 2 ( 3 ,oo) 12 8 V x/TT 


1 


+ o(t-^^ 


t > tQ,V > Vo ■ 


2V2 


3 / 1 , \ 3 , , , 

u+-lnt) <t< —-^{v + Int). 


272 
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Next, with fc = 1 in Lemma 4.6 and Proposition 4.7, 


lndet(/- 7i4rAi) = ^ ^ In |si| + Incg + In (l + tt ^2 q 

L^{s'^,oo) 12 8 V J \ J 

/3 1 


= ^ - In |si| + Inco + In (^1 + TT ^2 2 ^ 

We substitute this back into (4.15), 

lndet(/- 7 N:Ai) = ^ ^ In |s| + Inco + In fl + Tr-^-s- 

L^fs.oo) 12 8 V / 


+ ln(^l + 7r 2 2 + O {t 2^ 


and note that 


In I 


(^l + TT-h-Sti) +^V2(t-t'i) =ln(l + 7r-h-Jt-5ei^‘-’'^ + O , i\ < t < 4- 

Hence, 

Proposition 4.8. There exists to > 0 <^i^d Uq > 0 such that 


In det(/ — ^Kxi) 
uniformly for 


= :n5-xln|s|+lnco + lnn 1 + 77=2 


L2(s,C3o) 12 8 


t > to, V > Vo ■ 


j=0 




+ 0 




2 V 2 


(p + Int) <t< ( i; + - Int 


Remark 4.9. Note that the lower constraint on t is artificial: if we were to fix t < :^^iv + Int), the second 
factor in the product moves to the error term and we reproduce the result of Proposition f.l, after adjusting 
the error term. 

So far repeated integration has lead us to a sequence of results, 

(A.l) i-A tG [toiti] i-A tG [tljti] tG [£(^,£ 2 ], 

Prop. 4.3 Prop. 4.7 Prop. 4.8 

but this strategy can be continued indefinitely leading to the following result (here we also use Remark 4.9). 
Theorem 4.10. Given q G Z>i, there exists to = to{q) > 0 and vo = voiq) > 0 such that 


lndet(/ — yRTAi) 
uniformly for 


s" 1 
L2(s,Oo) 12 8 


9-1 


= W - ^ In |s| + Inco + In TT ( 1 + ^ 
I /> X n V VTT 

J=0 


O 


(H) 


t > to, V > Vo '■ 


3 3 

.v<t< —^ {v + q In t). 


2V2 ~ ~ 2V2 

5. Proof of Theorem 1.4 and asymptotics for eigenvalues 

The content of Theorem 1.4 is simply a combination of Theorems 4.10 and A.l: 

Corollary 5.1. Given x € K determine p G Z>o such that p = 0 for x < — 5 and X +5 <P + X+ i for 
X > — There exist positive to = to(x); I’d = 'i’o(x) such that 


D{JKi,l) = exp 


12 


p-i 


1=0 


f- 7^_| j_l 2^t_, 


4co]^(i+^2 2I It 1 ie3' 




uniformly for t > to,v > vo and Kj^. = f > fv^ — X^- la case p = 0, we take nj=o(- ■ ■) — 1 - 
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We now begin to derive asymptotic information for individual eigenvalues and employ techniques which 
have occurred previously in [7]. First, by positivity of e“'"Aj(l — 

det(/ - yiCAi) 


V € ^>0 • 1 “h G 


— V 

1 - Ap(s) 


> 1 = 


det(/ - yi^Ai) 


(s,oo) 


Now, we use Corollary 5.1 for the determinant in the denominator (with = x, x > 0 and thus p' = p € Z>o) 
as well as in the numerator (with x" = X + X ^ 0 and thus p" = p' + 1 = p + 1), i.e. 


> ( 1 




203 ' 




1 - Ap(s) 

uniformly for t > to,v > vq such that > \-\fi — X^j X ^ 0- Hence, after algebra. 


Vp G Z>o: 


Ap(s) 


> Cpt-P-sei^*, t^+oo; 


(5.1) 


1 - Ap(s) 

and the constant Cp > 0 can be chosen independent of v. Next, we make use of Lidskii’s Theorem: Valid 
for any i S Z>o, 

det(J - yi^Ai) 


det(/ - iXAi) 


(s,oo) 


= det(/ + e - K^s) 


n(r 

i=0 


A^(s) 


1 — A^(s) 


^ det (/ + e-"iX^(/-iX^)-i), 


with = i^Ai ■ Pi and projects on the space of eigenvectors of TXai with corresponding eigenvalues 
{Aj : j > P\. In this exact identity we can use the expansion of Corollary 5.1, i.e. as t ^ +oo,y f 1 such 
that > 1^2 - x^, 

p-1 


n(r 

j=0 


lL2-h-h-^-hAV2t- 


203' 


(l + O 


(5.2) 


£-1 

n(i 

i=0 




1 - Aj(s) 


^det (/ + e-"iX,(^-^^)”')- 


In here we first choose X = i-®- p = 2 and, say, £ = 1: 


VTT 


(i + o(H)) 


= I 


^-v ^o(s) 


det(l + e-"Ki(I-Xi)-^). (5. 


1 \ I ; V" ' /■ 0^-3) 

1 — Ao(s) / 

For §'/2t > V > — I Int we now see that det(J + e~^Ki{I — Ki)~^) > 1 together with (1.6) implies 

in (5.3) that 


^ie-iC 2 t 2 ^ 2 ^ lint 

^-^=0(1). i>i„. — . 

Next we choose x = | in (5.2), i.e. we take p = 3 and, say, 1 = 2: 


(5.4) 


1 H——2 4t 203^* ’ 
a/tt 




1+—2 2 it 

VTT 


= 1 + e 


-^ 0 ( 5 ) 


I - Ao(s) 


1 + e 


(i + o(H)) 

Jet (/ + e-”A',(7 - K^Y') (5.5) 


Multiplying through with the second summand of the first factor in the left hand side of (5.5), we obtain 


( /— 9 1 
1 + -v/tt 2Jt2 


g-§+2t+i 


-L _I_9 _3 2^/9^_. 

1H—=2 2 it 2e3v2t 

VTT 


(i + o(H)) 


(7 . e-.^) a. (7 . e-.7«7 - A,)-). 

















FROM GAP PROBABILITIES IN RANDOM MATRIX THEORY TO EIGENVALUE EXPANSIONS 


26 


For I's/Si — ^\nt > V > — | Ini, all factors in the left hand side are bounded, hence by positivity all 

factors in the right hand side have to be bounded and with det(J + e~'"— K 2 )~^) > 1, this leads us to 

2 /- lint 2 /- 31nt 

Tvy(;y = <^w; (>t„, 

After that, we let y = |, i.e. p = 4 and, say, £ = 3: After simplification in (5.2), 

T2Jt5e“s'^‘+"^ + v^25+Jt5e“5'^‘+''^ ^1 + ^ {l + O 

1 ^ z + . . _.^2t+v 7^2i + ltiAi(s) \ 

l-Ao(s) j l-Ai(s) ) 


= ^V7r2'it 

-y A2(s) 


X 1 + e 


det(/ + e-"iF3(/-£f3)-^) 


1 - A2(s), 

and by boundedness of the left hand side, also with the help of (5.1), for t>to,v> vq, 

TTy(:)=‘"w^ 5 ^-it^3^-20 

Iterating this approach for general x = P ~ = Z>i we derive a sequence of estimates. 


tt“*'2e 

1 - Aj (s) 


= 0(1), 0<j<p-2; 


2 e 3 * 

1 Ap_i(s) 


3A Int 


= 0(1); ,V2- p-^.y-^>>^.,>IV2- 


t 


l\ Int 


t 


Since these estimates are valid for any p G ^>i, we have in fact 

Aj(s) 


V j G Z>o : 


1 - Aj (s) 


<Djt ^ 2 03^*, i _|_oo 


(5.8) 


and the constant Dj > 0 can be chosen independent of v. With (5.1) and (5.8) thus 
Proposition 5.2. For any j G Z>o fixed, as t —>■ + 00 , 

1 - Ay(s) = + 0 ( 1 )), Cj > 0. 

We now recall a few general facts about the trace class operator Kai ■ ((s, 00 ); dA) O. From the identity 

Ai(A)Ai^(/i) — Ai(p)Ai'(A) 


ATai (A, /i) — 


A — /r 


Ai(A + t)Ai(t + fj,)dt 


it follows that Kai is positive definite with finite operator norm HlFAiH < 1 and trace norm 

pOO 

llAlAilli = / KAi(A, A)dA < ct, c> 0, t = \s\^. 

J S 

Note that for any positive definite, trace class operator B, 

I det(/ + B)-l\ < ||B||i exp (||.B||i), 

and thus, for any p G Z>i 

I det(/ + e-'^Kp{I - Kp)-^) - l| < e-'’\\Kp{I - Kp)-%exp {e-^\\Kp{I - Kp)-^\\,) 
But with (5.8), 

e-^\\Kp{I - Kp)-% < e-^\\KpU{I - Kp)-^ < e-"!^ < Cpti-^e 
With this back to (5.2) where k^. > , and we now let £ = p, 

det (/ + e-'’Kp{I - Kp)-^) =1 + 0 = 1 + o(l). 


(5.9) 

(5.10) 


VPo§ V2t-^ 


Summarizing, 
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Proposition 5.3. Given % S K, determine p G Z>o as in Corollary 5.1. There exist positive to = to(x)) '^o 
Vo(x) such that 


p-i 


1=0 


n 1+ 




p-1 


=n 1+ 


1=0 


l-Aj(s) 


(l +C> 


uniformly for t > to,v > Vq and > Iv^ — X^- here we take 11^=0 (■ ■ •) — 1 j (f P = 0 - 

At this point we use Proposition 5.2, take ln(...) of both sides and compare powers in t, this gives 


Cj = ^2-y-i, jeZ>o, 

VTT 


and Corollary 1.7 follows. 


6. Nonlinear steepest descent analysis associated with AT^eL “ part 1 
The Bessel kernel (1.4) is of type (1.10) with 

= Jai'/z), -ipiz) = JBess = (0,s), S>0 (6.1) 

and we will derive an asymptotic solution of RHP 1.12 for sufficiently large (positive) s and 7 close to 1 such 
that 

= ^ = 2 - 2 u = - ln(l - 7 ) > 0, t = s5; y e M>o, a G K>_i. (6.2) 


6.1. Preliminary transformations. We introduce the unimodular function 


4 'a(C) = V^e” 


4((e-‘"C)^) 

C)^^((e-‘X)^ 




-iiX4(e--C)^ 

-i(e--C)5iX^(e-- 




gifa^3^ CeC\[0,(X)) 


in terms of the modified Bessel functions I,y{z), K,j{z), cf. [33]. As before, we choose principal branches for 
: arg^ G (— 71 , 71 ]. Next we assemble 


argCG(0,f) 

vI/(C;a) = I, argC G (f, f) 

[(e-?), argCG(f,2.) 


(6.3) 


which leads us to the problem below. 

Riemann-Hilbert Problem 6.1. The parametrix '!'(('; a) defined in (6.3) has the following properties: 

( 1 ) 4'(C;a) is analytic for f G j U {0}) with 

fi = e*t ( 0 ,oo), f 2 = ( 0 ,oo), f 3 =e‘^( 0 ,oo) 

and all three rays are oriented as shown in Figure 5. 

(2) We observe the following jumps, 


4/+(C;a) = vI'_(C;a)(^_L;), C € Pi; vi.^(C;a) = vI'_(C;a)(/,„ J), C G Tg; 
^+{C,a) = 4/_(C;a)(_°,J), C G f2. 

(3) Near C = 0, in case a ^ Z, 


ik(C;a) = $(C;a)(e-'"C)®"^ 


2 sin 7ra | 

1 J 





argC G (0, f) 
argC G (f,f) ; 
argC G (f ,27r) 
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and for a G Z, 

vI/(C;a) = $(C;a)(e--C)^"^ (j 



argCG ( 0 ,f) 
argC G (f,f) 
argCG (f, 2 T) 


Here, ili((^;a) is analytic at (^ = 0 and we choose principal branches for fractional exponents and 
logarithms. 

(4) > oo, valid in a full neighborhood of infinity off the jump contours, 


iIi(C;a) = (e--C) 

X exp (e 


72 Vl 


U -Me- 


1 + 


8 (e-i’"C)® 


-(l + 4«7 2i \ 

2 i 1 -M 4a2 ) + 7 


where C,°‘ is defined and analytic for C, G C\(—c»,0] such that > 0 for ^ > 0. 



Figure 5. The oriented 
jump contours for the 
Bessel parametrix 4'(^;a) 
in the complex (j-plane. 



Figure 6 . “Undressing” 
of RHP 1.12 with s > 0. 
Jump contours of X(z) as 
solid lines. 


Remark 6.2. lUe observe that, cf. [33], 

= Jaiifh C G C\[0,oo) 

and hence for ( > 0, we have 

('I'aii(C))+ = Tre-*? J,(^/C), ('I^a2i(C))+ = v^e-*?^/CJ'(^/C)■ 

This enables us to rewrite the Bessel kernel in terms of the entries o/'I'a(C)j 

^Bess(A, M) = ^-, A, ^ e (0, S) 

in which all limits are taken coming from the upper half-plane. 

In our next move we reduce RHP 1.12 with (6.1) to a problem with ^-independent jumps and this step is 
the analogue of (2.5), (2.6) for the Bessel kernel determinant. Let s > 0 and set with help of Figure 6, 

{ I, z G Hi U H 2 U H 3 

(^_L°), ZGH4 
-GH. 

We arrive at the problem below: 

Riemann-Hilbert Problem 6.3. Determine X{z) = X{z; s,^; a) G such that 
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= = s + e'^( 0 ,oo), = (s,+oo), = s + ( 0 , oo) 

(2) The jump conditions are as follows, 


4 p(s) 
i=i ^ j 

C^) _ „ I p(^) _ I p(^) 


(1) X(z) is analytic for z G U {0,1}) with 


X+(z)=X_(z)(‘=-‘'“i-:), zef[^^-, X+{z)=X_iz){^,l), zGT 


is). 
3 > 


and 


x+{z) = x_{z){J^.°), x+(z)=x_(z)(^L;), 


(s) 


(3) Near z = s, with = {z € C : 3z ^ 0}, 


X{z) = X{z) [l+^{ ) ln(z - s) 


MM zgM+ \ i )’ ^ 

(- 10 ), arg(z- 3 )e(,r,f) 


I, Z G 


T 


else 


where X (z) is analytic at z = s and we choose the principal branch for the logarithm. On the other 
hand, near z = 0, 

1 — 


X(z) = X(z)(e--z)i'^ 3 (ii^)^ 


Xiz) = X(z)(e--z)i '^3 /_ (i_^)in(e--z)(oi) 


a g1>. 


(4) z —>■ oo, 

X{z) = (e-‘^z)-4‘"= 
with 




1 

^ VI 1 

1 fl i 


^ {/ + Xo,(e-*-z)-5 + O (z-i)} exp [(e-'-z) 5 a 3 


Xoo = -:r 


2 Vi -1 




lf-(l+4a^) 2i 


8 V 2i 1 + 4a^ / ■ 


6.2. Differential identity. The required identity is derived along the same lines as for D{JAi',j), cf. [15]. 
We skip details and summarize the final formula in the Proposition below. 


Proposition 6.4. For fixed 7 < 1, we have 

|lnD(Je..,7) = 0 = i 

in terms of the solution X(z) to RHP 6.3 and the limit is carried out for arg(z — s) G (f , 7 r). 


(6.4) 


7. Nonlinear steepest descent analysis associated with - part 2 
7.1. Initial transformation. We choose s > 0 and define 

T{z) = X{sz), zGC\(EtU{0,1}), 

where the contour Ep consists of four line segment, 

4 

Et = Pj.T; ri_7' = ( 0 ,l), r2,T = 1 + e ‘3 ( 0 ,00), Pa^r = (l,+oo), r4_T = 1 + e*^( 0 ,00); 
i=i 

oriented “from left to right”, see Figure 6 . This leads us to 

Riemann-Hilbert Problem 7.1. Determine a function T{z) = T{z;s,^;a) G which is uniquely 

determined by the following properties: 

(1) T{z) is analytic for z G C\(St U {0,1}). 
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(2) We have the jump conditions 

r+(z)=T_(z)('^-;"“p), zefi,T; T+iz)=T_iz){\l), z G Pg.r; 

and 

r+(z) = r_(z)(^_L.°), zGr2.T; t+(z) = r_(z)(^.L °), zGr4,T. 

(3) The singular behavior near z = 0 and z = 1 is unchanged from the one stated in RHP 6.3, modulo 
the change of variables T{z) = X{sz). 

(4) As z ^ oo, with t = 




^ e-‘?'"3{/ + Xoo(e-‘’^sz)-5 +C>(z-i) jexp t{e-^^z)^a3 


7.2. Normalization transformation. We introduce the ^-function, for z G C\K, 


g{z) = -(e-‘’^(z- 1))® +yin 


l + (e-‘’^(z-l))5 


V = 


^1 - (e-‘^(z - 1))2 J t 

where C“ as well as InC are defined and analytic for (j G C\(—oo, 0]. We also choose 

'l + (e-*’"(z-l))3' 


(7.1) 


—TT < arg 


< TT. 


1 — (e *’^( 2 : — 1 ))2 

The relevant analytical properties of the g-function are summarized below. 


Proposition 7.2. The function g{z) defined in (7.1) is analytic for z G C\((—oo,0) U (l,+oo)), more 
precisely, 

g±iz) = -Vl - + y In ^ ^ 


and 


9± 


\J\ — z — 1 
1 + v^l — z 


=F iiry, 2 G (— 00 , 0); 




G (0,1); 


as well as 


Also for z —)■ 00 , 2 ^ K, 


, 1 — Vl ~ zj ’ 

g±iz) = ±iV 2 - 1 ± iParg , 2 G( 1 ,+oo). 

(e-‘’^2)5 + g{z) = iTrP + (^V - 0 (e-*’^2)-5 + O ( 2 "^) . 

In the next transformation we set 

S'( 2 ) = e-‘’^*'''^ 3 ^( 2 )e‘®(^)'^^ 2 G C\(Et U {0,1}) 

and are lead to the problem below. 

Riemann-Hilbert Problem 7.3. Determine S{z) = 5 ( 2 ; 5 , 7 ; a) G such that 

(1) 5 '( 2 ) is analytic for 2 G C\(TtT U (— 00 ) U {0,1}). 

(2) The limiting values S±{z),z G are related by the equations 


0 


^ g—i7rag2i(;(2) 1 / ’ ^ ^ , *^+('^) *^—(^) y Qi7Ta^2tg{z) ly ’ ^ ^ ^ 4 ^ 7 ^, 

and 


1 


5+( 2 ) = S'-( 2 ) 
as well as 


0 1 

-1 0 


2 G Ts^t; S+{z) = S_(2)e-2"*'"'^^ 2 G (- 00 ,0); 


S+(2) = S_(2) 


,-i7ra g-t(j<B 0 as+ 23 ( 2 )) 

0 e'^“ 


, 2 G Fi^t- 
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(3) Near z = 1 and z = 0 with S{z) analytic at either point, 

eiTTtVaa = S{z) [-^ + ^ () ln(z - 1) 


qI 2 *^*^3 


{\l), zgh+ 
I, z G El” 


arg(z-l)G (f,T) 
arg(z - 1) G (tt, ^) , z^l, 
/, else 


and, as z ^ 0, z ^ Et, 


S'(z)(e ‘’"sz) 2‘^3 ( Q 2sin7roj^ a^Z 


iirtV(T3 


S{z)e 


-tg(z)a3 _ 


= ^(z)(e-‘"sz)5‘^3 


^(l-7)Me *’'sz)([ji) 


, a G Z. 


(4) As 


S{z) = (e-‘’^sz)-J'^^e-'’^‘^'"^^ {/ + S^{e-^^z)-^ + O (z-^)} 

with 

+t(^2V- a^. 

As before, the following observations are crucial: With 0 < r < | fixed, 

3?(g(z)) <0, z G (f 2 ,T U r 4 ,T)\^(l, r) 
and for z G (0, l)\£)(0,r) and sufficiently large t > to,v > vq subject to (6.2), 


>fBes= + 2ff (^) = “ 2^1-^ + ‘2V lu 


>5>Q. 


(7.2) 


(7.3) 


Hence, we expect the major contribution to the asymptotic solution of RHP 7.3 to arise from the line 
segments (— 00 , 0) U (0, 1) U (1, + 00 ) and two small neighborhoods of z = 0 as well as z = 1. 

Remark 7.4. We will again represent tV = x as 

X = fc + a; fc G Z>o, 

and hence the jump on (—oo,0) in RHP 7.3 equals 

.5+(z) = 5_(z)e-2"‘“"^ z<0. 

7.3. Analysis of model Riemann-Hilbert problems. The outer parametrix, for z G 


P(“)(z) = (e-‘’^s(z - 1)) 

with the Szegd function 

l + (e-‘’^(z-l))' 


,±(1 -A Ae-^-(^-l))^+l ^ 

•\/2 yl 1/ I (e“*’^(z — 1))5 — 1 


(I?(z))"^ (7.4) 


I?(z) = 


z G 


2q: 


x>(z) = e*’"“ 1 + -—:—^ +0{z ^) , z^oo; 

(e“*’^z)2 


1 - (e-i’^(z - !))■ 
displays the properties summarized below. 

Riemann-Hilbert Problem 7.5. The parametrix P^°°\z) has the following analytical properties 
(1) PA)(Z) is analytic for z G C\M and we orient the real axis from left to right. 
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(z) = , z G (0,1); 


pr\^) = pr\z){ °^h), zG(i,+oo) 


(2) The limiting values P^\z),z G M are square integrable and related by the jump conditions 

(3) As z ^ oo, z ^ M, 

For the local parametrix near z = 1 we use ideas from [ 6 ], Section 3.4. Consider the unimodular function 


S(C) = 2 VT 


1 ^ A o\ / (cA 


, CgC\(-oo,0] 


(7.5) 


vO iy \^-cAi7A)'(cA cA^r)'(cA, 

which uses the Hankel functions and and the principal branch for A : arg^ G (—7r,7r]. 

Riemann-Hilbert Problem 7.6. The function B{Q G defined in (7.5) has the following properties: 

(1) i?(C) is analytic for ( G C\(—oo,0]. 

(2) Provided we orient the negative half-ray from —oo to the origin, we have 

s+(C) = i3-(C)(J " 2 ^), Ce(-oo, 0 ). 

(3) 7lsC^0,C^ (-00,0], 


1 


5(C) = oV^(o ilAC) 


1 0 


InC f-1 -1 


27ri \ 1 1 


with the principal branch for the logarithm and 

' Jo(cA(l-fln2) _j„(^i)(i+|ln2) \ 

-c^J^(CA(1 - f In 2 ) - fJo(cA J^(cA(i + f In 2 ) - fJo(cAy 

-00 (-K)'' 


B(C) = 


2 i 


--E^A+ 1 ) _ 


fe =0 


1 


-2k -2k) {k\y 


1 \ (- 1 )'^ 


ICI < r; 


'^o(C=)— Z)o (fc !)2 

A^^(cA = Er 2 fc^ 


(4) The function B{(j) is normalized so that as C ^ 00 , ^ ^ M, 


B{C) = C 




72 VI 1 


8 C^ 


- 2 i -1 


+ 


1 4i 


128C V“4i 1 


O 


(rs)}e‘<- 


(7.6) 


In terms of B{(j) we can now define the local parametrix as follows: consider 


gi,rtVa3p(i)(2)e-t3C)<T3 ^ E^^\z)B{C{z)) 


I- 


1 — 7 /—1 —1 


27ri 


1 


1 )ln(C(^)) 


1 2 RCTs 


(°M z G 11+ [(e-"“i)’ arg(z l)G(3,7r) 

arg(z-l)G(7r,f); 0<|z-l|<i 

else 


Z ^ . 


(7.7) 


where we make use of the locally conformal change of variables 


Ciz) = 7e- g{z) - - In 


l + (e-‘’^(z-l))5 


t V ^ ~ (e“‘’^(z — 1))2 




1 + 


4fc(z — 1) 
3(t- 2fc) 


+ 0((z-l)2)}, 


so that 


, X _i2L^ 3 I / N Ct ^ / 1 + (e ‘""(z - 1)) 2 

(j2(^z)=e ‘ 2 isgn(! 3 z) < g(z) --In ' 


t \ 1 — (e“''^(z — 1 )) 2 




3 5 
4’ 4 


(7.8) 
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and the locally analytic multiplier, for 1^ — 1| < 

f I 0 - 


to -1 g JJ+ 1 

(^1 0 j, ^ g-lfa<T3glfo-3_ 


1 1 


my 


I, 2 G H" V-1 ly 

It is easy to verify the properties of P^^\z) which are summarized below: 

Riemann-Hilbert Problem 7.7. The parametrix P^^^z) defined in (7.7) has the following properties: 

(1) is analytic for z G D{1, U {I}). 

(2) Since the change of variables ( >—>■ ((z),z G P(l, |) is locally conformal, we obtain from RHP 7.6 
directly the jump behavior 


py\z) = 


(1)^ 


1 


0 


py\z) = p!i’{z) 


(1)^ 


g-i7rag2tg(2) ^ 

1 0 
gi7rag2tg(2) ^ 


, 2: G r2,T 17 £* ( 1, - ) , 


1 


, 2 : G r 4 _j’ n D 1 , 


and 


P^Hz) = J) 


( 1 ,- 1 ; 


py\z) = py\z) 


,-ma g-t(j<Boas+2g(2)) 

0 e'^“ 


, ^ e ( -,i 


This matches exactly the jump behavior near z = 1 of RHP 1. 3. 
(3) Near 2 = 1, with the principal branch for the logarithm, 


gi7rtVa3p(l)('^)g-ig(2)<T3 




0 2^ 


(yi), zeU+nD{i,l) 

Z gU- nD{i,l) 


I, 


(^-La°), arg( 2 :-1) G (f,7r) 

(_j,a°), arg(0 - 1) G (tt, ^) , 2^1, 2^ Et, 

^ I, else 

and we have thus the same singular structure near z = 1 as in RHP 1. 3. 

(4) In case t —)■ +oo,7t 1 subject to (6.2), we obtain from (7.6) and (7.8), 


pC>(.) = A~)(4b+^ML)^_2ie- 


1 


-2ie-'™ 

-1 


(M(2)) 


-r 


128C(2) 


(7.9) 


X M{z) 


1 

_4iei7ra 


4ie 


-1 


j(M( 2 )) ' + 0(t-3)|; Miz) = {V{z)) 


ZGH+ 


z ^ . 


uniformly for 0<ri<|2: — 1| < r 2 < Here we use that on the annulus 


exp 


± 2ick2) =0{t°°), t-)>+oo, 7tl: = 2 - 2 [x + ^ 


Inf 

t 


Remark 7.8. Observe that from properties (2) and (3) in RHP 7.7 we obtain 

5(2) = iVi(z)pW(z), 0<|2-l|<i, 

with Ni{z) analytic at z = 1. 

For the remaining parametrix near z = 0 we shall use again classical orthogonal polynomials: define for 
k G Z>o, with principal branches for fractional exponents, 

L((;a)=( f. y I (e-‘’^C)®'"^e-5‘^'"^ CgC\[0,oo), (7.10) 

V7fe-r9/c-i(C) Jo Qk-iPJt ^ t=c/ 
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in terms of monic Laguerre polynomials = 0 = 7 fe-i, cf. [33]: 

poo 

^(C) = C^'+ Ofc,fc-iC^ ^ + a/c_fc_ 2 C^ ^ + O ^) , C —>■ oo; / q^j\t)q^j,\t)t°'e *dt = hkSjk, 

Jo 

with 

k 

hk = k\r{l + k + a), keZ>o', ak,k-i = —k{k + a), ak,k -2 = ^{k + a){k — l){k + a — 1), fc G Z>o. 

Remark 7.9. The normalization in [33] of the Laguerre polynomials {L^^\c.)}kei->o cigain different from 
the one we choose: we have to use the relation ql^\C) = (— C £ 

Riemann-Hilbert Problem 7.10. For any k G Z>o,a > —1, the function L{(f;a) G defined in (7.10) 
has the following properties 

(1) L{f; a) is analytic for f G C\[0, +oo) and we orient the positive half-ray from 0 to +oo. 

(2) Along this ray we have 

' 3 —iTTa ^ 


L+{f\a) = L_{C;a) 


0 


C G ( 0 , 00 ). 


(3) Near C = 0, 


L{C-a) = L(C;a)(e-‘"C) = "^ (o 2 sin™J, ^ ^ Z; 

L(C;a) = L(C;a)(e-'X)t-3 
with L{f; a) analytic at ( = 0. 




a G Z; 


(4) C ^ oo,C ^ [0,oo), with jk = - 7 ^ 

1 




_ 2 

^2 I _ 


^ 7^ ) +0 (C”') \ 

Ik 2 2 7 fc+i J \ 


(7.11) 


The required model function near z = 0 is now given by 


where 


P(0)(z) = ii;(°)(z)L(C(z);a)e5-Be..-aeis(.).3^ ^ e P (^0, J) \ , (7.12) 

p(o)(z) =p(~)(z)(l?(z))-"^(e-‘"C(^))-5"^tt-3(^(^))-3^ 0< |z| < i 


with the choice 


m = ciz) 


l + (e-‘’"(z-l))5 


fx = 22 ^=^“ 


1 — -z + O (z^) , z —)• 0 , 


1 — (e“‘’^(z — 1 )) 2 

is analytic at z = 0. In (7.12) we have employed the locally conformal change of variables 


CL) = 2, U)-bn' ' + 


r) + 1 ] — :n|l + -+ (2 ( 2 ^) 1, 71 < 


t \ 1 ~ (e“'^(z — 1 )) 2 

and we collect the relevant properties of P^°^ (z) below. 

Riemann-Hilbert Problem 7.11. The parametrix P^^\z) has the following analytical properties 

(1) P(°)(z) is analytic for z G P(0, 7 )\(- 7 , 7 ). 

(2) With local analyticity of ( = C,{z), 

p|o)(z) = P®(z)e-2"“^'^3^ 


p|°)(z) = p 1 °)(z) 


( -^, 0 ) ; 


g-iira g-t(><BeBs+2g(z)) 


z G 0 
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which matches the jump behavior of S{z) in a neighborhood of z = 0, compare RHP 7.3. 
(3) Near z = 0, compare RHP 7.10, 


^i7rtV(T3 


^ r 'n:a 

gi^tva3p(0)(^)g-is(z)a3 ^ P® (z)^ In (e-‘"sz) (g J) 


a e Z; 


i.e. we also have a matching of the singular behavior with the one stated in RHP 7.3. 

(4) As t ^ +c »,7 f 1 subject to (6.2), we derive from (7.11) (using also that M>o 3 x = ^ + ct)) 

( 7 fc'r^(^)(e-'-t-ic(^))-“' 


p(0)(z)=p(~)(z)<{/+—(p(z)) 


1 


7fc_i/32(z)(e ^'^t 


7k-1 
7k 


{'D{z)y 


e{z) 


(P{z)) 


-3 ( i(:l.p-\z){e---t-K{z))-^\ 

2l.vMo-i'n’y.-l/‘I.r'Oa 1 7fc-l I \^\Z)) 

2 7fe+l / 


-7fe_2/32(z)(e '’"t ^C(z))‘ 


O 


(t-3+2|“l) 


(7.13) 


uniformly for 0 < ri < |z| < r 2 < 7 . 

Remark 7.12. From properties (2) and (3) in RHP 7.11 we deduce 


S{z)=N^{z)P^°\z), 0<|z|< 


where Nq{z) is analytic at z = 0 . 


We have now completed the necessary local analysis and can move on to the comparison of the explicit 
model functions P*'“^(z), P*'^^( 2 :) and P^^\z) to the unknown S{z) of RHP 7.3. 

7.4. Ratio transformation and first small norm estimate. Use (7.4), (7.7), (7.12) and define 

{ (p(o)(z))“^, |z| < r 

(pW(z))-\ |z-l|<r ; uj = -s^N{S^-{a + 2 a)a 3 )N-y^^ 

(P(°°)(z))”\ |z| > r, |z-1| > r 

with 0 < r < 7 fixed and we put, compare RHP 7.3, 


_ -nraas 


N = e 


NO -iU-‘ 
1 


e“‘4'^3gi’racr3^ 



Figure 7. The oriented jump contours for the ratio function R{z) in the complex z-plane. 


Riemann-Hilbert Problem 7.13. Determine R{z) = R{z; s,'y; a) G such that 
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(1) R{z) is analytic for z S CySn with square integrahle boundary values on the contour 

Efl = dD{Q, r) U dD(l, r) U (r, 1 — r) U ^(r 2 _T U r 4 _ 7 ’) n{zGC: |z—1|> r}^ 
which is depicted in Figure 7. 

(2) We have jumps R+{z) = R-{z)Gii{z), z G Sjj with Gr^z) = Gfi{z] s,j; a) and 

Gr(z) = z G dD{0,r); Gr{z) = pW(z)(p(“)(z))-\ z G 80(1,r); 

followed by 

Gfl(z) =p(-)(z) {P^°^Hz))-\ z G f2,T n {z G C : |z - 1 | > r}, 

Gfl(z)=p(-)(z)(^^.^,^i^(,) (p(-)(z))-\ zGf4.Tn{zGC: |z-l|>r}; 


as well as 


GR[z) = pr\z) 


(oo). . (l , I. 


0 


1 


(Pi°°^(z)) , zG(r, 1-r). 


Note that by construction, see Remarks 7.8 and 7.12, there are no jumps inside D{0,r)U D{l,r) and 
on (—oo, —r) U (1 + r, +oo). Moreover R{z) is bounded at z = 0 and z = 1. 

(3) z —>■ oo, 


P(z) = 


1 0 
w 1 


-ivr ',-io'3 -iirtVcrs 


(e ‘’^sz) 


1 






V2 Vi 1 

X (p(“)(z))”^ =/ + G((e-‘’^z)-5) . 

Observe that we obtain from (7.2), (7.3) and (7.9) the following small norm estimate for 

GR(z;s, 7 ;a) = s^’^^GRiz; s,^-a)s~i'^^, z G Sr. 


(7.14) 


Proposition 7.14. Given x G R> 0 )a > ~1 there exist positive to = toiXT^)No = vo{x,a) and c = c(x,a) 
such that 


||GR(-;s,7;a) - 7||L2nL”(SH\9D(o.r)) < -, Vt > to, '9 = -ln(l - 7 ) > no : r = 2t - 2 (^x + 2 j 

For the remaining contour dD{0,r) we have to be more careful: From (7.13), uniformly for z G dD(0,r), 
Gr{z-,s,x] a) - I = (e"‘’"(z - i^Rf {z)t~^~''^°‘ + Rf {z)t~'^~'^‘^ + R 2 {z)t~^ 

+ P+(z)t-2+2“ + p-( 0 )t- 2 - 2 “ + R^{z)r'^ + o (t-3+2|“l) }(e-‘’^(z - 1))J'"3 

with t-independent coefficients 

Rti^) = i7fe (] e™; Rf{z) = -i^,_i^^(z)e-™ 


2C(z) 


2C(z) 


1 1 
-1 -1 


and 

2C2(z) 

all of rank one as well as 
RNz) = 

C(^) 


1 

2C2(z) 


1 -1 


1 1 
-1 -1 


7. R4{z) = 


4C2(z) 


/ 7fc-2 _|_ 7fc-i 7fc-2 _ 7fc-i \ 

iTratTs I 7fc 7fc+i 7fe 7/e+i J ^i-TratJa 

I 7fc-2 _ 7fc-i Ik-2 I 7fc-i I 

V 7fc 7fc+i 7fe 7fe+i / 


We have introduced /l(z) = t “/3(z), ^( 2 :) = iC(2^) and 

'1 + (e-‘^(z-1))5 


7(z) = 


1 — (e“'’^(z — 1)) 2 


rC(^) — 2^“< 1 —-z + G (z^) >, z —>■ 0. 
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Similar to the approach for the Airy kernel we use also here matrix factorizations, 

Gij(z;s,7;a) = + R2{z)t~^ + (R^iz) - R'l{z)R2{z))t~'^^'^°‘ 

+ R^{z)t-^-^‘^ + {R^{z) - R+{z)R^{z))t-^ + O (t-3+2|“l) }(e-‘-(z - 0 < a < ^; 

Gr{z- s,7;a) = E-{z){e-^^{z - + R+{z)t-^+^^ + + {R^{z) - i?r 

+ i?+(z)t-"+2“ + (i?4(^) - R^{z)Rt{z))t-^ + O (t-3+2|“l) }(e-‘-(z - -^ < a < 0 

with the invertible meromorphic factors 

E+{z) = (e-‘"(z-l))-^"^(/ + i1■+(z)^-l+2“)(e-'"(z-l))3'^^ 

E-{z) = (e-‘"(z-l))-J"3(/ + i?r(z)f^""“)(e-*"(z-l))J'^b 

We have thus 

GR{z-s,r,a) = z S dD{Q,r) 

and can move on to the next transformation. 


7.5. Singular Riemann-Hilbert problem and iterative solution. Put 


Q{z) = R{z) < 


-i‘"3if+(z)s5'^3. 

14 

< r 

5 

14 

> r 

-i'^^E-{z)si'^F 

14 

< r 

5 

14 

> r 


, 0 < a < 2 


-5 < a < 0 


(7.15) 


and obtain 

Riemann-Hilbert Problem 7.15. Determine Q{z) = Q{z; s,'y; a) G such that 

(1) Q{z) is analytic for z G C\(IJr U {0}) with square integrable boundary values on the contour T,r 
depicted in Figure 7. 

(2) The jump conditions are as follows, 

Q+iz) = Q-{z)GR{z-,s,j;a), z S r); 

g+(z) = Q_(z) i<“<t 577(0,r). 

-5<a<0 

(3) We have a first order pole singularity at z = 0, 

^1 


Q{z) = Q{z) 


\z\ < r 


with Q(z) analytic at z = 0 and 

l2-4fc-2a-r^-r+2Q 


a+= 


1 — i7^ l2“4fe-2a-2^-l+2. 

as well as 
T+ = 


a G 


2a 


0^ X ; <^ = 


-i7fe_i24'=+2“-it-i- 

1 + i7fc_i24fe+2a-2i-l-2a ■ 


a G 




e27riagi Q 






T- = 


1 1 

_e277iagi Q 


, OL G 


4'“ 


(4) As z ^ 00 , 


Q(z) =/ + o((e-‘-z)-!) . 


In order to complete the nonlinear steepest descent analysis we perform a final transformation, define for 
z S C\Si?, 

g(z) = {z/ + R±}L(z)i (7.16) 

and obtain the problem below 
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Riemann-Hilbert Problem 7 . 16 . Determine L{z) = L{z', s,j; a) G such that 

(1) L{z) is analytic for z G with the contour shown in Figure 7. 

(2) We have identical jumps as in the latest Q-RFIP 7.15, i.e. 

L+{z) = L_{z)GQ{z;s,'y;a), zGEr. 

(3) L{z) is analytic at z = 0 provided that 

= a±L(o)r± (^55 (r±)-'|L(o)-a±L'(o)r± J) ■ 

(4) As z —>■ oo, we have L{z) —>■ I. 

At this point we use a similar factorization as in (7.14), 

GQ{z-,s,-c,a) = zGT^r 

and summarize our results. 

Proposition 7 . 17 . Given x = A: + a S K>o with k G Z>o,—^ < a < ^ as well as a > —1, there exist 
positive to = Ao(X)«)) 'f'o = i'o(X) a) o,nd c = c(x, a) such that 

\\GQ{-,s,r,a) - /||L 2 nL~(aD(o.r)) < Vi > to, v = -ln(l - 7 ) > wo : u = 2t- 2 (^x+ Int. 

Hence, from general theory [21], the latter Proposition together with Proposition 7.14 implies 

Proposition 7 . 18 . Given x G ffi>o and a > —1 there exist positive to = to(x, o),iio = vo{x,a) and c = 
c(x,a) such that RHP 7.16 is uniquely solvable for t > to,v > vq : p = 2t — 2 (x + §) Int. The solution 
L = L{z; s,j;a) can be computed iteratively via the integral equation 

1 /* d 

L(z) = /+— / L_{w){Gq{w) - I) — zgC\T.r-, L(2;s,7;a) = s3'^3L(2;;s,7;a)s"3'^^ 

2711 Jsh ' w-z 

using the estimate 

l|i-(-;s,7;a) - yt>to,v>vo: u = 2t - 2 (^x + ^) Int. 


8. Extraction of large gap asymptotics 
We employ the central differential identity (6.4), i.e. 


d 


F(JBess;7) = ^lni7(JBess;7) = - ^6*"“ (A"^ (z)X'(z)) . 
and the limit is taken with arg(z — s) G (f ,7r). 

8.1. Asymptotics for the derivative. Tracing back the transformation sequence 

X(z) I—)■ T(z) I—)■ S(z) I—)■ R(z) I—)■ Q(z) i-A L(z) 
we derive for (6.4) the following exact identity 


E(JBess;7) = - 




IP—>-1 


76 


27ris 


{(P(i)(«;)) '(P«(ip))'} 


21 


27ris 

^giTra 

w->-i 27 ri 5 


(S Xn')S'(w))^ 


10—>-1 


|(P^^^('u;)) ^R ^(w)R'(w)P^^X'a’)J 


It;—vl 


The first summand is derived from RHP 7.7 


ye 
27ris 


7ra ^ 

^{(P(i)(u;))- (P«(rc))'}^ 


,1 a k k, ,, 
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and for the second we have from (7.15) and (7.16), 


7e 


27ris 


|(P^^^('u;)) ^(r(;)i?'(i(;)P^^^('u;)| 


7e 


fi 27ris 




7e 


27ris 




= fl _ ~-2 

2s V t 


^ e-2-“|L-i(l)B±L(l)}^^ + £ “ t ) 

Note that with Proposition 7.18, for z S 

Hz) = I+^(f {GQiw)-l)^^ + ^ (f (GQ(ii;)-/)^^+o(t-2-4|“l) 

2711 JdD(0,r) ^ ' w- Z 2-K\ JoD(l,r) ^ ' W- Z \ ) 


ldD(fi,r) 

and hence by standard residue computations. 


L(l) = |'^-mi„{ 2 - 2 |a|. 2 }^ |gi,ra^3^ 

L'{1) = e-'’^“‘"^ + O min{ 2 - 2 |a|. 2 }^ |gi 7 ra<T 3 


where 


Aa 




8 V t 


4a 


1 2fc 


7l2 = 1 - — 


-1 


t ) V-4a2-^ * 


3t-2k 


+ 


l 7fc- 

1 ( 

70 

4) 

' 7fc 

1 

V4 

oj 

7/c-i / 

'0 

l^ 

1 

7fc \ 

^1 

oy 

1 


and 


ih _ r)=t4fc —l± 2 a f f 4^4^ I ITA;— 1 ; (4~) . tzt _ /tih 


^ 0=2 


A± — _ .4± 

zn1 1 / 1 • -1 / I ’ ^2 — ^0 ■ 

^4 -4/ 17;, , (-) 


We now obtain 

{l-1(1)L'(4)}^^ =e2’"*“s5{A±t-l-2|“l +242f^+C>(f™"42-2|a|.2}^ 


= e^^ss 


a 


21 

' /fc. 


-+ -(A: + a)±2=^'*"-^=""“t 

2t ' 


±4fc-l±2a.f-l-2|a| J ITfc-l) ( + ) 


and with 


I + O 


+ O f 4 -™“{ 2 - 2 |q|, 2 }'\ 


40-3 i7ra(T3 


5 4 ^ e 


also 


= ^±g-i7ra<T3g-i 


L±1(0) = 

i<^3 |/ + c> {/ + C> I sJ 


jt73gl7rQ;fT3 


= e2-‘“s5a±|(y^ T4 


so that 

[l-\1)B^L{1)} 

We summarize (using that 7 = 1 — = 1 + 0 {t~°°)), 

. . 1 a (T^ k k 

(JBess;7) = + + 

and can now move to the integration. 


(/fc + a) + o(t- 2 - 2 |“l) 
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B.2. Integration of expansion (8.1). Start with k = 0 and 0 < a < |, i.e. 


^ 4 2v^ 4s 2yi ^ l+r(l + a)7r-i2-2«-3t-i+2a 

and make use of 

Lemma 8.1. Let t = t{s) = ss and 0 < Sq < s. For any f G L^(so, s), we have 


( 8 . 2 ) 


Since 

we have that 

and thus in ( 8 . 2 ), 
1 

~2 


= - = 1 - 


/(t(u))du = 2 / wf{w)dw; to = Sq ^ 

Jto 

2-2 u =-ln(l - 7 ) > 0 , x = fc + i 


j2oL j^—2k — a^—v-\-2t 


2do(a)i 


e 


fa^t{u))^= f 

J Sq J Iq 

= In +i(l + a)/ - 

llJ — tn tn 


1 + do(a)ru“^““e“"e2“ 
t 


2 do(a)w^”“e-’'e2“' dw 
+ do{a)w~^~°-e~^e'^'^ w 


with doia) = r(l + a) 7 r“^ 2 “^““^. Suppose t > to,v > vq are sufficiently large such that 

0 < 2 (I' + ol’^t) < t < d- (l + a)lnt), 

and thus the base point of integration equals Sq = j(t + alnt)^ or equivalently to = 5(11 + a In t) >0. Then 

1 


■ [ CT'''(t(w)) ^ = In (1 + do(a)t ^ “e^* ”) + O (t 
Jsn M2 V / 


and we can integrate ( 8 . 2 ) as follows, 


In 


det(/ yit^Bess) 1 ^ 2 ( 0 ^^) 
^det(/ - yi^Bess) 1^2(0, so). 


1 


= -t(s - So) + a(yi - V^) --r In 


4 

+0 (t-i 


+ ln(l + do(a)t-^-V‘-") 
(8.3) 


where we use that 


f C>(t- 2 - 2 “(w))dM= / C> dw = / C> (w-i+“e- 2 ’"e'') dw = C> (t-i) 

«/ Sq ^ tQ J t'Q 

uniformly as t —>■ + 00 , y t 1 such that 0 < ^(m + alnt) < t < \{v + (1 + a) Int). But with (B.l), we know 
that 


lndet(/ - y if Bess) 


1 


G(l + a) 




since > 2 — a^. Substituting the latter back into (8.3) we have thus 

Proposition 8.2. For any given a > —1 there exist tg^a) > 0 and vo(a) > 0 such that 


lndet(i - y if Bess) 
uniformly for 


s ^ n , 1 , 

= --+ ay's -- In s + In Ta + In 

L 2 (o,s) 4 4 


(^1 + do{a)t 


— l — a^2t—v 


o{t 


- 1 ') 


i ^ V > vq ■ 0 < - (v + a In t) < i < - (^J + (1 + a) In t) 


with do (a) = r{la) 7T ^2 


— 1 o —2a—3 
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In the general case k G Z>o we let 


and first derive the following two expansions. 

Lemma 8.3. For ik <t<i'f.,i.e. a G [—^,0], 

/det(/-7ifBess)L.(0,,) \ 1 ~ ^ , I /- 

VdetU-7^Bess) 1^2(0,,,)/ 4 4 


+ 2fc(v^ — ySk 


- 2^k + a)lii 


+ In (l + i7fc_i24'=+2“-2t-i+«+2fce-2‘+^^ - In (l + i7fc_i2^'=+2“-2) + O (f^) . 

Proof. We use that 

I [ a-{t{u))^ =lnfl + i7fc_i24'=+2“-2i-i+“+2fce-2‘+"') - In f 1 + i7fe_i2‘‘'=+2“-2') + O (f^) 

and the stated identity follows now from (8.1). 

Lemma 8.4. For i'^.<t< ik+i, *-e. a. G [0, ^], 


□ 


In 


det(/ 7^Bess) I ^ 2 (Q g 


- 2^k + a)\n 
Proof. Note that 


det(/- 7i4:Bess) 

s 


—) =-J(s-Sfc) + afv^--\/^) Y + 2fc('v^- 

,si)/ ^ 2 k \ / 


L2(0. = , 

^2^ _ j..^-l2-4fe-2a-2j.-l-a-2fe„2t-j) 




o{t 


-11 


cr+(t(M))^ = In ('1-17^12 
^ Js' U2 V 


— 1 o—4/c —2 a—2^— 1 —a—2fcg2i —? 




and substitute this back into (8.1). 


□ 


At this point we follow the logic carried out for the Airy kernel determinant, i.e. first choose fc = 1 in 
Lemma 8.3 and use Proposition 8.2, 


In det(J — yi^Ai) 


L2(0,; 


Si 

' 4 


+ ay/Il - ^ Insi + InTa + In (^1 + do(a)^ + O (t 


so that 


lndet(/— yATeess) = —- + aVs — — Ins + InPo + In fl + do(n)i ^ “ 

L^(o,s) 4 4 V 

and we used that 


— l — a^2t — v 


2(t-ti) = -ln(ti+“e-2‘+^) . 

Summarizing, 

Proposition 8.5. Given a > —1 there exist tQ^a) > 0 and Vo(a) > 0 such that 


lndet(/- yATBess) 

uniformly for 


L2(0,s) 


=-^ + ai/s-^ Ins + InTa + In ^1 + do(a)t ^ “ 


— l — a^2t — v 


- (?; + (1 + a) Inf) < t < -{v + {2 + a) Inf). 


O 


Inf 

t 


- 01 ^^ 


t > to, V > Vo ■ 
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Now we continue with k = 1\ use Proposition 8.5, 


L2(0.si) 

let di{a) = r(2 + and thus back in Lemma 8.4, 

s 

L2(o,s) 4 


lndet(/ - ■jKBess) + a^/§[ - ^ Ins'^ + InPa + In (^1 + do{a)t^ + O > 

ack in Lemma 8.4, 

2 

+ a^/s—— Ins + InTa + 2(t — {[) + In ^1 + do{a)t^ 


In det(/—7^86 
+ ln| 




But since 


2(t-t;) = ln 

we obtain in fact 

Proposition 8.6. Given a > —1, there exist to{a) > 0 and vo(a) > 0 such that 

^ + a^/s -^ In s + in Tq -f 

j=o 


In det(/ - 'yKBess) 
uniformly for 


L2(o,s) 


2 ^ / 
ly/s — ^ Ins + liiTfl + In -\- O i 

')=n \ 


“) 


t > to, V > Vo '■ - (i; + (2 + a) Int) < t < ^ + (3 + o) ■ 


Remark 8.7. Also here the lower constraint on t is artificial since for t < ^{v + 2\nt), the second factor of 
the product only contributes to the error term and we reproduce Proposition 8.5 with adjusted error estimate. 

At this point we would iterate the procedure (just as in the case of the Airy kernel) and obtain in the end 
the following result 

Theorem 8.8. Given q G Z>i and a > —1, there exist to = to{q,a) > 0 and vo = vo(q,a) > 0 such that 

2 9-1 


2 H~ ^ / 

lndet(J - ^Kbbss) ^ ~ X ^ IT ^ 


x) 


with dj(a) = j\T[l + aj)T: ^2 uniformly for 

t>to, V > Vo ■ 0 < -(p + alnt) < t < - (p + ( 2(7 + a) Int). 

9. Proof of Theorem 1.8 and asymptotics for eigenvalues 

To obtain Theorem 1.8 we combine Theorem 8.8 with expansion (B.l), 

Corollary 9.1. Given y € K and fixed a > —1 determine p G Z>o such that p = 0 for x < o,nd 
X+l <P<X+i for X > — 5 . There exist positive to = io(X) o),vo = wo(X) o) such that 

-D(JBess; 7 ) = exp -^+a^/s s~^°‘^ Ta\Y[l + dj{a)t~'^^~'^~°-e^*~^'^ ^1 + O ^max ^ 

with dj{a) = j\T{l + a-\-j)TT~^2~^^~'^°‘~^, uniformly for t > to, v > vo aird = y > 2 — 2(x + f )^. In 
case p = 0, we take n 7 =o(' ■ •) = 1- 

At this point the logic is just as before we only have to use that is positive definite with finite 

operator norm ||Arg“gg|| < 1 and trace norm 


l41lli= r<ess(A,A)dA<ct, Vt>to, t = 

Jo 
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Appendix A. Refining Corollary 1.2 in [7] 

In this appendix we will improve the error estimate derived in [7], Corollary 1.2.. In more detail, we shall 
prove the following Theorem. 


Theorem A.l. As s —>■ —oo, 7 t 


det(/ - 
uniformly for 


(s,oo) 


= exp 


S'" 

12 


|s| sco(^I + C>(^s Co = exp 


-In2 + C'(-1) 


(A.l) 


ln(l — 7 ) . 2 
(-s) 


^ 

2 O 


Expansion (A.l) was proven in [7] with a weaker error term. The result in loc. cit followed from the 
Tracy-Widom representation of Z?( Jai; 7 ) in terms of a Painleve transcendent and the derivation of transition 
asymptotics for the latter. 

Here, we will apply the method of integrable integral operators and use RHP 1.12. In more detail, (A.l) 
will follow from a comparison of the nonlinear steepest descent analysis of RHP 1.12 for 7 = l,s —)■ —00 
(corresponding to = + 00 ), cf. [15], to the nonlinear steepest descent analysis of the same RHP for 
7 t l,s —)■ —00 such that is fixed. We will see that the difference between xr^. > fixed and 

xr^i = +00 is “almost beyond all orders”. 


A.l. Comparison of nonlinear steepest descent analysis. We go back to RHP 1.12 tailored to (2.1) 
and observe its validity for any choice of 7 < 1. The same also applies to the first two transformations 

Y(z) hA X(z) HA T{z) 

carried out in Section 2.1. After that, we can formally copy the g-function transformation (3.1), however for 
the scale 


-^ 2 , +00 


fixed, 


or 


XT,, = +00 


we shall see that setting H = 0 in (3.1) is sufficient. With this, the next transformation 

5'o(^) =T(z)e*(3(^)l^=«)'"^ zeC\YT 

leads us to the following problem. 


Riemann-Hilbert Problem A.2. The normalized function Sq{z) = So{z;s,j) G is characterized by 

the following properties 

(1) So{z) is analytic for z G C\(Et U {0}). 

(2) The limiting values S±{z),z G are related by the equations 

5o+(z) = 5o_(^)(^^ 2,,(,V^„ J), ze (r2Ur4)\{0}; So+iz) = So-iz) J), zePaUO}; 

and, along the entire positive half-ray, 

/l p-t(>^Ai+23(z)+=o)\ 

-S'oa( 2) ='S'o-(z) ( Q ^ j, z€(0,+oo). 

(3) The singular behavior is identical to the one stated in RHP 3.4 subject to the replacement V = 0. 

(4) Also the asymptotic normalization at z = 00 E-r is the same as in RHP 3. 4, modulo V = 0. 

Note that, as before, with 0 < r < | fixed, 

^giz)\v=o) <0, z G (r 2 U T4)\D{0, r) (A.2) 

and, since 

g{z)\v=o > ^ S (O’ + 00 ); 9 lv=o = 
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we have now for x G (0,+oo)\D(0,r), 

+ 2g(z)\v=o > 0 with equality iff 



and 



(A.3) 


as a consequence of our choice > | \pi. But this means that we only require a simplified local Riemann- 
Hilbert analysis near z = \ when working with the scale G [f's/S, +oo]. In fact, for ^Ai G [|v^,+oo) we 
take A: = 0 in (3.12), i.e. 







CG 


and adjusting all subsequent steps of the construction to V = 0 = Q; = fcis sufficient. The matching (3.16) 
is then replaced by 

P(i)(,) = |j+ ^(0 1) + o (i-f ) I . + x) ^ 

uniformly for 0 < ri < ji: — ^| < r 2 < On the other hand, for = +(X) we do not need any separate 
analysis near z = \. 


Remark A.3. If we were to keep x’^j > fixed, there would be no need for a local analysis near z = ^ 
at all, see (A.3). The difference between > I's/S and xr^j = +(X) is thus “beyond all orders” and between 
^ ^Ai — “almost beyond all orders”. 


The local model functions in (3.5) and (3.10) can be used again, subject to the replacements V = 0 = 
a = k. 

A.2. Small norm theorem. In the given situation the ratio transformation consists of the replacement 

V -1 I 




V—O—OL—k'’ 


z < r 


V=0=a=k' 2 !^^ 


w = -|s|5(A5ooA^-^), 


lv=o=a=fc> \z\>r,\z-^\>r 
when xr^j G [Iv^,+c») is fixed and without P^^\z) 


,, „ , for X,. = + 00 . In either case 0 < r < | and 

V—Q—(y.=k o 


we use 


N = 


1 




(' be-'L 

V2I-1 V ■ 

We are lead to the problem below which is formulated on the contour shown in Figure 4 for G [fv^, + 00 ) 
and in Figure 8 for xr^j = +c». 



Figure 8. The oriented jump contours for the ratio function Po(-z) in the complex z-plane, 
case xr^j = + 00 . 


Riemann-Hilbert Problem A.4. Determine Ro{z) = Ro{z;t, G such that 
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(1) Ro{z) is analytic for z G with square integrable boundary values on the contours Srq, 

dDio, r) u ((r 2 u r 4 ) n {z e C : |z| > r}) ; 

Sro = dD{0,r) U dD ^ + t, oo^ U ((Fa Ur 4 ) n{z e C : |z| > r}); 

shown in Figures 4 <§• 

(2) On the contour we have jumps Ro+{z) = i?o-(z)GRj,(z; s, 7 ), z G Er with 
Gr„(z) =p(0)(z)|^^p(p(~)(z)|^^J-\ z e dD{0,r); 

Gr„(z) ^ e dD Q,r) 


Gr„(z)=p(-)(z) 

Gr„(z)=p(-)(z) 


I g-t(>JAi+2g(z)|v=o) 

lv=olo 1 


(P(“)(z)|^^j \ zG (r, --r ) U (r +-,+00 


v=o (e^Jlv^o J) -e(r 2 Ur 4 )n{zGC: |z|>r}. 


These jumps hold for >c^. G + 00 ). For xr^. = +00 we simply have to dispose of dD{^,r) U 

(r, i - r) U (r + 5 , + 00 ). 

(3) z —>■ 00 , we have Ro{z) —>■ /. 

It is evident that the latter problem admits direct asymptotical analysis, in fact from Proposition 3.15 
and (A. 4) we obtain at once 

Proposition A.5. There exist tg o,nd c positive such that 

I|GRo(-G, X^Ai)-^llL2nL“(ER„) < -T> Vf > to, >fAi > liGRo(’;i,+Oo)-^l|L2nL~(EHo) ^ IT’ '^^>^ 0 - 

C 6 O 13 

This estimate implies unique solvability of the Pq-R-HP in A^(Srj,). Moreover, the solution satisfies 


Ro{z]t,>c^fj = |s| 


t 2 

TTR 


|4'^^ t^oo,7tl: >fAi > 3^ (A-S) 


Po(z;t,+oo) = |s| + ^ t^oo 

uniformly for z G C\Er(,. We shall now consider 

x{z;t,Kj^J = Ro{z;t,>iAi){Ro{z;t,+oo)) \ z G C\Er; t>to, G 

Because of (A.5) and (A. 6 ) the function x(z;t, xr^J is well defined and it again satisfies a RHP 

Riemann-Hilbert Problem A.6. Determine x(z) = x(,z:;t, G such that 

(1) x{z) is analytic for z G C\E^ with square integrable boundary values on the jump contour 


(A. 6 ) 


•\/ 2,+00 I fixed. 


1 


E^ = (9P(0, r)U5P -,r U r, - -r U r + +oo 


1 


1 


shown in Figure 9. 



Figure 9. The oriented jump contours for the ratio function x{P} in the complex z-plane. 
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(2) The limiting values are related via the jump conditions X+(^) = x) with 


Gx{z)=Ra-{z-,t,+oo) P^°\z)\ {P^°\z)\ y=o ) Rol{z;t,+oo), zedD{0,r); 


>c—-\-oo 

- 1 ' 


Gxiz) =Ro{z;t,+oo) P^^\z)\y^^{P^°°'>{z)\y^^^) ^ i?o ^(z; t,+oo), zGdD(-,r 


and for z G {r, ^ — r) U {r + ^, + 00 ), 


Gxiz) = Ro{z;t,+oo) 


p(°o)(^) 


'1 Q-t(^Ai+^9G)\v=,o) 

lv=olo 1 


(p(-)(z)|^^o) 


Rq +oo)- 


(3) As z ^ 00 , we have x(z) —>■ I- 

Using (3.10) together with (A. 6 ) we obtain the following important small-norm estimate 


Proposition A.7. There exists to and c positive such that 

16 6 

But this implies again unique solvability of the y-KHP for t > to, x > and we have 

= + ^eC\E^, 00 . 

We summarize the relation between the nonlinear steepest descent analysis with fixed S [fv^, + 00 ) and 
= -|-oo in the following Corollary. 


Corollary A.8 . There exists to > 0 o,nd c > 0 such that 

t / 1 \ N 

t 2 


Ro{z;t,xJ = \s\-z-^ ->1 + 0 


1 + kl 


s\*'^^Roiz;t,+oo), yt>to, x^Ai G 


-a/ 2,-boo) (A.7) 


uniformly for z G C\E^. 

A. 3. Proof of expansion A.l. We recall the differential identity (2.8), 

^lnP(JAi;7) = ^(X-i(z)W(z))^^[^^ 

in which the derivative is taken for any fixed 7 < 1. Repeating the steps carried out in Section 4.1, here 
however tailored to the transformation sequence 


X(z) !-)■ T(z) H> So(z) H> Ro(z), 
we derive immediately the following exact identity 

F(JAi;7) = ^lnP(JAi;7) = ^ 7 //- |{i?o-'(0)i?U0)}^^. 

We can now either directly use the Pq-R-HP to compute the remaining matrix entry or simply refer to 
our previous computations of Section 4.1. If we choose the latter we have to keep in mind the additional 
transformation Q(z) 1 —>■ L(z) and tailor it to the current situation. In either case, we obtain after indefinite 
integration with respect to s, 

lni:)(JAi; 7 ) = ^ ln|s|+lndo+ C> -boo, 7 f 1 : G 

where do > 0 may still depend on 7 (but not s) at this point. 

Remark A.9. Observe that for > I's/S all power like error terms in (A. 8 ) are in fact x-independent. 
This matches the result of Remark 7.4. in [7]. 


72, +00 j (A. 8 ) 
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It is clear that 


rfo( 7 ) = Co (1 + 0 ( 1 )), 7tl; Co = exp 


-^ + C'(-i) 


but we would like to have a better error estimate. To this end, note that 


^ In £)( Jai; 7 ) = -tr ((/ - "/KaO ^Kai) 


1 


(s,oo) 


poo 

/ R{X,X\t,>c^,)dX 

J S 


with, (cf. [15], Section 2), 




, Ae(s,oo); (') = —. 
^-^-A dz 


But from (A. 7), we get, after tracing back the transformations, 

where we used the classical gap asymptotics from, say, [15] in the last equality. Hence, back in (A. 8 ), we 
deduce 

^\ndo{^) = O (t-iy 

and thus 


ln(io( 7 ) = Inco + O t->+oo, yfl: 

which completes the proof of (A.l). 




^ a / 2 ,+00^ , 


Appendix B. Refining Theorem 4 in [20] 


This part of the appendix complements the previous one and we derive the analogue of Theorem A.l for 
the Bessel kernel. In more detail, we shall prove 


Theorem B.l. For any fixed a > —1, as s ^ + 00,7 f 1, 


det(/ - 7 A:Bess) 
uniformly for 


L2(0,s) 


= exp 


= -- 


ln(l - 7 ) 


•+( 

1 + 0 ( 

> 2 - 

Int 

a -. 


t 

■^Bess 

= +00 


C?(l + a) 

( 277)5 


(B.l) 


Another proof of the gap expansion (i.e. = + 00 ) appeared in [25] using operator theoretic methods 

valid for jSftaj < 1. Here we will use once more the method of integrable integral operators and RHP 1.12 
subject to ( 6 . 1 ). 


B.l. Comparison of nonlinear steepest descent analysis. Since the initial RHP 1.12 is valid for any 
7 < 1 we can again apply the first two transformations 

Y{z) HA X(z) HA T{z) 

as worked out in Sections 6.1 and 7.1. After that we set V = 0 in (7.1) and define 

So{z) = z e C\St. 

This leads us to the problem below 

Riemann-Hilbert Problem B.2. The normalized function So{z) = So{z;s,^) G is characterized by 

the following properties 

(1) So{z) is analytic for z € C\(St U {0,1}). 
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(2) The jump conditions are as follows, 


*5^0+(■2-) i7rag2i(;(2)|v'^o l) ’ ^ ^ ^2,T: ‘^0+('2') Sq—(^z) ( i-7rag2i(;(2)|l) ’ ^ ^ ^4:,T 


and on the positive real axis, 


Sq+[z) = So-{z) 


,-iTTa g-t(>CBesB+2s(z)|v=o; 
0 6*’"“ 


, ZGTi^t; So+{z) = So_{z) 


0 1 

-1 0 




(3) The singular behavior is identical to the one stated in RHP 7.3, subject to the replacement V = 0. 

(4) The normalization at z = oo is identical to the one stated in RHP 7.3, also here subject to the 
replacement V = 0. 

We have, as before, with 0 < r < | fixed, 

^g{z)\v=o) <0, zG (r2,TUr4,T)\^(l,r), 

but since 

5'(2)|v=o >-1, 2:G(0, 1); g(0)|y=o = -l, 

we now also have that for given a > —1 there exist to > 0, 0 < r < j such that for z G (0,1)\-D(0,r) and 
t > to, 

+ 2g(-z)|v=o > 2 — 2^/1 — z — > S > 0. 

This means we only need a simplified local analysis near the origin z = 0. In fact, take fc = 0 in (7.10) 

J_ +ap-t_^N 
^ 271-i Jo ^ i-C 


L(C;a) 


k=0 


0 


1 


-C J ^ g C\[0,oo), 


and adjust (7.12) to 17 = 0 = a = fc. The same adjustments have to be made for the local model functions 
in (7.4) and (7.7). 


B.2. Small norm theorem. We are left with the final transformation: put 
1 0^ 


Rq{z) = 

with 0 < r < i and 


r (p(°)(^)) 

, ,^5o(z) (P«(z))-^ 

[(p(“)(z)) 


Y—0—ct—k ’ 

1^1 

< r 

Y—0—ct—k ’ 


— 1 < r 

1 Y—0—ct—k ’ 


> r, \z 

1 



N = ^ 

i 

B e“' 

V2 

VI 

ly 


cu = 


-s^^ {N{S^ - aa3)N 

(B.2) 


The jump conditions of Rq{z) are posed on the contour shown in Figure 7 when ^rBess + ^ [2,oo). For 

■^bobs = have no jumps on dD{0, r) U (r, 1 — r). In either case we easily derive the following estimate 

for 

GRoiz;s,-f;a) = si’^^Guoiz', s,-f; a)s~i'^^, z G Sro 
with Grj,(z) denoting the jump matrix associated with (B.2). 


Proposition B.3. Given a > —1 there exist to = to(a),Vo = vo(a) and c = c(a) such that 

c 

l|GRo(-;'S>7;a)--^IU 2 nL~(EflQ) < p Vt > to, i'= - ln(l - 7 ) >'^o : v>2t-ah\t. 

From this estimate we obtain unique solvability of the Pq-R-HP in L^(I]rq) and the solution satisfies 

a > -1 : s"^'"^Ro{z)s~^'"^ =1 + 0 , t-J>oo, yfl: x^bbbb > 2 - (B.3) 

uniformly for z G C\T,r„. The required comparison follows from writing Rq = Ro{z;t, >fBeBBi analyzing 

= Roiz;t,K^^^^-,a){Ro{z;t,+oo;a)) \ z G C\Er(,; t > to, x^bbbb + G [2,+ 00 ) fixed. 
Riemann-Hilbert Problem B.4. Determine x(z) = x(z; t, o) G such that 
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(1) xi^) *■5 analytic for z € with square integrable boundary values on the jump contour 


= dD{0, r) U dD{l, r) U (r, 1 — r) 


shown in Figure 10. 



Figure 10. The oriented jump contours for the ratio function Ro{z) in the complex 2 ;-plane. 


(2) We have jump conditions X+i^) = niith 

Gxiz) =Roiz;t,+oo;a) ^ Rg^(z;t,+oo;a), zG 

Gxiz) =Ro-{z;t,+oo;a)\p^^'>{z)\.^ {P^^\z)\ y=o )”^] +oo; a) 


=+°° 


and for z G (r, 1 — r), 
Gxiz) = Ro{z;t,+oo;a) P^°°\z) 


I e“‘™e“‘(^Bes.+2s(^)lv'=o: 

v"=o V 0 1 




(3) Near z = oo, the function behaves as x(z) —)■ I. 

Now we use (B.3) and derive for the conjugated jump matrix 


dD{0,r); 
z G dD{l,r); 

\ ^ Pq ^(z; t,+oo; a). 


the estimate 

Proposition B.5. For any given a > — 1 there exist positive tg = to(a) and c = c(a) such that 

-^IU^nL“(s^) < p Vt>to, > 2-ai^. 

This provides in turn the analogue of (A.7) for the Bessel kernel and the remaining argument is just as 
in Section A.3. 
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